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Abstract 

We study the collective excitations of holographic quantum liquids formed in the low 
energy theory living at the intersection of two sets of D-branes. The corresponding 
field theory dual is a supersymmetric Yang-Mills theory with massless matter hyper- 
multiplets in the fundamental representation of the gauge group which generically live 
on a defect of the unflavored theory. Working in the quenched (probe) approximation, 
we focus on determining the universal properties of these systems. We analyze their 
thermodynamics, the speed of first sound, the diffnsion constant, and the speed of 
zero sound. We study the influence of temperature, chemical potential, and magnetic 
field on these quantities, as well as on the corresponding collisionless/hydrodynamic 
crossover. We also generalize the alternative quantization for all conformally AdS^ 
cases and study the anyonic correlators. 
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1 Introduction 


The understanding of new phases of matter is one of the main goals of fundamental physics. 
At low temperature and non-zero density, Fermi liquids are commonly described by the Lan¬ 
dau’s phenomenological theory [^, in which the elementary excitations take place near the 
Fermi surface and are fermionic quasiparticles. Landau’s theory has been very successful in 
dealing with an ample variety of low temperature materials [^. However, there are examples 
of systems whose behaviors are not well described by the Landau theory. The gauge/gravity 
duality is a new principle which could shed light on in establishing new paradigms in 
the case of systems with strong interactions and without a quasiparticle description. In 
particular, there is some hope that holography could help to classify compressible states of 
matter, he., states with non-zero charge density which varies continuously with the chemical 
potential p. 

In this paper we adopt a top-down approach to this problem and explore the properties 
of matter engineered with intersections of D-branes of different dimensionalities at non¬ 
zero density and low temperature. We will consider generic intersections in which N^. Dp- 
branes intersect Nf Dg-branes, with q > p, along n common spatial directions. We will 
denote by (n |p T g) such a D-brane intersection (for example (3 | 3 T 7) for the familiar 
supersymmetric D3-D7 conhguration). The gauge theory dual of this (n | p T g) conhguration 
corresponds to a (p-l-l)-dimensional SU{Nc) gauge theory in which one adds Nj fundamental 
hypermultiplets living on a (n -|- l)-dimensional defect [^. 

Our description will be valid in the’t Hooft large-W limit, with large ’t Hooft coupling. 
In addition, in our approach we will assume that iV/ W and we will treat the Dg-branes 
as probes in the gravitational background created by the Dp-branes, which corresponds to 
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the quenched approximation on the gauge theory side. 

By turning on a suitable worldvolume gauge held we will add a non-zero baryonic charge 
density . In addition we will switch on a magnetic held along two of the spatial directions 
of the worldvolume. The embedding of the brane probe is characterized by a function, 
which measures the distance between the two sets of branes. The gauge theory dual of this 
distance is the mass of the hypermultiplet helds. In the present paper we will only consider 
massless fundamentals for which the embedding function is trivial. We leave the analysis of 
the massive case for a future work. 


We aim to determine universal properties for this holographic matter, which do not 
depend much on the particular intersections and are common to all of them. It turns out 
that the different observables studied depend only on the dimensionality p of the bulk theory, 
as well as on an index A, dehned in terms of p, q, and n as: 

A = 2n + ^ (p — 3) (p -I- g — 2?7, — 8) . (1.1) 

This index A is the same for several intersections (for example, A = 6 for all (p | p T (p -|- 4)) 
conhgurations). Both p and A determine a universality class in the set of intersections 
analyzed. 


Among the observables we study are the thermodynamic properties (entropy and specihc 
heat), speed of hrst sound, as well as the excitation spectra. The latter can be obtained 
by looking at the quasinormal fluctuation modes of the system. At zero temperature we 
look at the holographic zero sound, which is a collective mode hrst found in and then 
generalized to hnite temperature in and to magnetic helds in [^. The holographic zero 
sound has been on focus in many diherent contexts 10-30 and, in particular, for some 


of the D-brane intersections included in our analysis. Our results reproduce the values 
found previously for the speed and the attenuation of the zero sound and generalize them 
for arbitrary intersections. The zero sound mode is the dominant one at sufficiently low 
temperature, where quantum ehects dominate over the thermal ehects. 


At enough high temperature thermal ehects become more important than quantum ehects 
and the system enters into a hydrodynamic regime. We will hnd that the dominant collective 
mode in this regime is a dihusion mode and we will be able to calculate an analytic expression 
for the corresponding dihusion constant. At temperatures between the collisionless and 
hydrodynamic regimes there is a crossover transition which we will study numerically for the 
diherent values of (p. A). 


One of the main objectives of this paper is to study the inhuence of the magnetic held 
on the collective excitations of the diherent systems; previous studies include D3-D7’ [^, 
D2-D8’ 23 , D3-D5 20 22 , and D3-D7 22 . It is known that, for sufficiently large B at 


T = 0 and non-zero density, a quantum phase transition takes place in the D3-D7 model 31 


Similar transitions may happen in other intersections, we thus assume that B is below the 
critical value of the possible phase transition]^ 

= 6 models: Sakai-Sugimoto, D3-D7’, and D2-D8’, are subject to modulated instabilities, so they 
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In the presence of the magnetic held B parity is broken and the longitndinal and trans¬ 
verse oscillations of the brane probe are coupled. This coupling corresponds to the mixing 
of the dual held theory operators that is induced when 5 7 ^ 0. For small values of B we 
will be able to decouple the equations of the zero sound mode. The result of this analysis is 
that the zero sound is now gapped, with the gap being given hj B/fi for any values of (p, A) 
at T = 0. This result generalizes the ones, e.g.^ in [9,22 
theorem 


and is in agreement with Kohn’s 


32 . At non-zero temperature there is a critical magnetic held above which the 


zero sound becomes massive [^. As in [^, the location of the collisionless/hydrodynamic 
crossover is insensitive to the strength of the magnetic held. We will also study the dihusion 
mode and determine the magnetic held dependence of the dihusion constant. 


When the intersection is (2-|-l)-dimensional there exists the possibility of imposing mixed 
Dirichlet-Neumann UV boundary conditions [^, he., to perform an alternative quantization, 
to the quasinormal modes. This alternative quantization corresponds to rendering the charge 
carrying excitations to anyons (particles of fractional statistics) and is characterized by some 
constant n, which measures the degree of mixing of the boundary conditions. We will study 
the collective excitations of the anyonic huids in the presence of the B held and we will hnd 
that the zero sound is generically gapped, although the parameter n can be hne-tuned to 
produce gapless spectrum precisely when the anyons feel a vanishing ehective magnetic held. 
This parallels the hndings 33 34 in the anyonic superhuid formed in the D3-D7’ model in the 
incompressible phase. We will also study the dihusion constant and the DC conductivities 


of the anyonic systems. Related works in the D3-D5 model appeared in 35,36 


The rest of this paper is organized as follows. In section we present our general setup 
and analyze its thermodynamical properties, both at zero and non-zero temperature. In 
section!^ we begin to study the huctuations of the probe. In section [^ we concentrate on the 
case of vanishing magnetic held. We hnd analytic equations for the dihusion constant and 
the spectrum of the zero sound. These analytic expressions are compared with numerical 
calculations. Section is devoted to the study of the ehects due to the magnetic held in the 
dihusion constant and the zero sound. In section [^ we perform the alternative quantization 
for the (2 -|- l)-dimensional intersections. Section [^ contains a discussion on the scaling 
behavior of diherent dimensionful quantities. For example, we argue that A determines the 
scaling dimension of the charge density. Section contains our conclusions and discusses 
future lines of research not included in the present work. 


The paper is completed with several appendices that complement the analysis made 
in the main text. In appendix we collect the equations of motion for the huctuations. 
In appendix we solve the indicial equation for the huctuations around the horizon. In 
appendix we hnd the correlator of two transverse currents, in the absence of the magnetic 
held. Finally, in appendix we detail the method used to solve the coupled equations of 
the R 7 ^ 0 zero sound. 

need to be studied at large enough temperatures (relative to charge densities). 
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2 Y)p-Y)q systems with charge and magnetic field at 
non-zero temperature 


In this section we introduce our setup and study some of its properties. Let us consider 
a generic ten-dimensional background metric at zero temperature (we will soon relax this 
assumption) of the type: 


*10 = 3u{r)dt^ + g„(r) + ... + (dx^f] + g„{r)dg■ dff , 


( 2 , 1 ) 


where y = (j/^,... are the coordinates transverse to the Dp-brane and the functions 

9tti Qxxi and depend on the transverse radial direction r (r^ = y-y)- We note, that these 
metric components have explicit functional forms in the case of (black) Dp-brane background, 
but we retain from revealing them until Section 2d^ to stress that the following treatment is 
rather general. We now embed a Dg-brane probe extended along the directions 




( 2 . 2 ) 


We will refer to this configuration as a (n | p ± g) intersection (n is the number of common 
spatial directions of the Dp- and Dg-branes). This intersection is represented by the array: 


X 

Dp : X 
Dq : X 


X X 
X — 


XP 

X — 

— X 




X 


y 


.9—p 


We will denote by z the coordinates y transverse to the Dg-brane: 


9+n-p-q 


) . 


(2.3) 


with z™' = y'^ for m = l,...,9-|-n — p — g. Moreover, let us define the coordinate p as: 

p2 = {y^f + ... + {y^-^f . 

Since: 

{dyf = dp^ + + dz 
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(2.4) 



(2.5) 

as: 




.. + {dxPf 


+ dz"^ . 


(2.6) 


We will consider embeddings with x^^^ = ■ ■ ■ = x^ = constant and \z\ = L, with L also 
constant. Then, = p^ + L? and the induced metric on the Dg-brane worldvolume takes 
the form: 

dSq+i = gtt{r)dt^ + g 3 , 3 ,(r) [{dx^f + ... + [dx'^f] + g^^(r) [dp^ + ] . (2.7) 
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In what follows we will consider massless embeddings with L = 0 (we will retnrn to L 7 ^ 0 
case elsewhere). In this case the r and p variables are eqnal. Moreover, we will switch on 
a non-zero temperatnre, which amounts to including a blackening factor fp, in such a way 
that the induced metric becomes: 

do^ 

dsj+i = gtt{p) fp{p)dt^+9xx{p) [{dx^f + ... + {dx^f]+grr{p)[jj^+p^dnl_^_^] . ( 2 . 8 ) 

Let us now compute the DBI action of the Dg-brane with a non-zero worldvolume gauge 
held F with components along p, t. Moreover, when the number n of Cartesian coordinates 
on the Dg-brane worldvolume is greater than or equal to 2, we will allow a non-zero constant 
magnetic held B along the directions x^ and Thus, we will take F to be given by: 

F = A[dpAdt -I- Bdx^ A dx^ , (2.9) 

where A[ = dpAt and we have chosen a gauge for A such that Ap = 0. The DBI action for 
this conhguration becomes: 

Sng = -Td, j ^/-dei{g + F) = j dp^^grr\gu\ - , (2.10) 

where lR(n,i) is the volume of the {n + l)-dimensional Minkowski space, (j) is the dilaton of 
the background, Af the constant: 


Af = TDgVo\{S'^-^-^) , 


( 2 . 11 ) 


and H is the function: 


H(p) ^ g-p (sL + B") . (2.12) 

The gauge held component At is a cyclic variable in the action ( |2.10[ ) and has an associated 
constant of motion. Then, we can write: 


d 


Vha', 

\/drr \gtt\ 


where d is a constant (the charge density). Inverting this relation we get: 


(2.13) 


_ d \J Prr \gtt I 

‘ “ ~7WTW ■ 

By using this result we obtain the on-shell Lagrangian density: 


^{on—shell) rr Idtt I jj^ 

Cns, 


(2.14) 


(2.15) 
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2.1 Intersections in the Djg-brane background 


Let us now write some of the equations for the background corresponding to a stack of 
Dp-branes. The metric and dilaton in this case are given by: 


9tt 9xx 




= I - 


(7-p)(p-3) 

2 


where i? is a constant radius. The blackening factor fp is given by: 


fp = l- 


er ■ 


where Vh is the horizon radius, related to the temperature T as follows: 


T 


7 — p 
dvr 


5-p 

/y^ 2 


(2.16) 


(2.17) 


(2.18) 


In the following we will scale out the constant R or, equivalently, we will directly take 
R = 1. In this case, the function H corresponding to a (n|p T g) intersection takes the 
form: 


H ^ p + p'*”-'' B‘‘ , 


(2.19) 


where the constant exponent A is the combination of p, g, and n written in (1.1). We 
will show below that the numbers p and A determine the thermodynamics and collective 
excitations of the intersection. 


2.2 Thermodynamics at zero temperature 


We begin by studying the thermodynamics of the intersections at zero temperature in the 
absence of a magnetic field by following 37 . According to the standard AdS/CFT dictionary. 


the chemical potential p at zero temperature for the brane intersection is given by the 
boundary value of At(cx)) and for Dg-brane probes entering the Poincare horizon can be 
written as an integral: 


/i = At(cx)) 




y/ 9rr \9tt I 

+ (P 


dp . 


( 2 . 20 ) 


When the metric is given by (2.16) and H = (see (2.19)), the chemical potential becomes: 


where 7 is a constant 


p = d 


dp 






= , 




r(i + 



( 2 . 21 ) 


( 2 . 22 ) 













The on-shell action of the probe (without the Minkowski volume factor f^{„,i)) is given by 
the integral of the on-shell Lagrangian density (2.15): 


Nort¬ 


on—shell 


= -Af 


P 




dp , 


(2.23) 


which is divergent and must be regulated. We will do it by subtracting the same integral 
with d = 0. We get: 


qreg = - Kf 

^on—shell 


A 


A 

p2 


p^ + d"^ 


dp . 


(2.24) 


The grand potential density = fl(/i) is given by minus the regulated on-shell action: 


o — _ 

“ ~ ‘^on-shell ' 


After performing explicitly the integral in (2.24) we arrive at: 

2 


n = 


A + 2 

From we can obtain the density p as: 


Af d^~^^ = 


A + 2 


Af J 2 p +2 


dfl , 

Moreover, the energy density e is given by e = Q + p p: 

A 


e = 


A + 2 

The pressure is just p = —and therefore: 

2 


Af 'y d^^ ^ . 




(2.25) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


We notice that A is thus related to the polytropic index of the equation of state. Finally, the 
speed of (hrst) sound is dehned as: 


o dp 2 

“•"S = A- 


(2.30) 


which follows immediately from the relation between p and e. We wish to warn the reader 
though, that obtaining the speed of sound at T = 0 is a little bit of a stretch, as being 
beyond the applicability of the hydrodynamics. 

We have thus found that the zero temperature thermodynamic behavior is determined 
uniquely by the index A dehned in (1.1). In subsection 2.5 we will discuss which models 
violate the bound on the speed of sound as expected for strongly coupled held theories with 
a gravity dual. 
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2.3 Thermodynamics at non-zero temperature 


We now study the thermodynamics of the probe at T 7 ^ 0 and S = 0. The chemical potential 
is given by: 


h = 


dpA[ = d 


'rh 


'rh 




11 1 

dp = 1 + -; 


2’ A’ 


A’ d2 


The grand potential at T 7 ^ 0 is given by: 


= -S:T.UeU=l^ 


A 


"Th 


X 

P 2 


xj p^ + 


- 1 


dp 


2 




A + 2 


(A + l)d 


F -,i + t;2 + t;--^ + 


A 


A’ d? 


2M 

A + 2 


(2.31) 


(2.32) 


The last term in (|2.32|) is independent of the density and therefore of chemical potential. 

= (2.33) 


We dehne the p-dependent part Af2(p) of the grand potential as: 

2M 


A + 2 


In order to study the behavior of the system with the temperature, let us consider the 
case in which is small. Expanding the chemical potential (2.31) at next-to-leading order 
in Vh, we get: 


A+i 


The expansion of Af2 is: 

AQ = 


+ = 7ci--r. + ^^3^^ + 0(r, 

■—^— Af'yd^'^^ -—— 

A + 2 ^ (A + l)d 


2 A+l^ 


Plugging the expression of d (|2.34|) into (|2.35|) we End: 
2 


An = 


A + 2 


2 (^p + rh 


,A+l^ 1+A 


2(A + 1) d2 


Af 


(A + l)d 


^A-hl 

h 


The entropy depending on the density is given by the following derivative: 


s = 


d{An) 


dT 


d{An) 


From (2.36), we obtain 


dvh 


d{An) 


dvh 

dT 


5 — p \7 — p 


Ett 


5-p ^^2=3 c}(Ar 2 ) 

_/ 5-p 


dvh 


dvh 


"-vd-Ap + .... 


(2.34) 


(2.35) 


(2.36) 


(2.37) 


(2.38) 
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Therefore, we get: 


s = 


An (^ . 1 (A'kT\^ 


5 — p \7 — p 


M dT^-p 1 + 


2d^ 


W-pJ 


+ 


Let us now compute the specihc heat cy from the formula: 

^ dT 


(2.39) 


(2.40) 


For p 7 ^ 3 we only need to keep the leading term in (2.39) to hnd the behavior at low T, 
which only depends on p. We hnd: 

2(p — 3) / An \ ^ 


cv{p ^ 3) = 


UdT'^v + 


(2.41) 


(5-p)2 V7-p. 

Notice that cy is linear in T (as for the Landau-Fermi liquid) only for p = 4. For p = 3 the 
entropy is instead: 

s(p = 3) = nUd + ^ ^ + ■ ■ • , (2.42) 

and the specihc heat depends on A in the following form: 


Cy{p = 3) = + 


(2.43) 


Notice that A = 2n for p = 3, as follows from (1.1). Notice also that the entropy at T = 0 
is non-vanishing in the p = 3 case, as pointed out in [^[^. This degeneracy suggests that 
there is an instability towards a non-degenerate ground state, see, e.g., 38 . However, to 


date such instabilities have not been found. As discussed in 39 , the backreaction of the 
havor D-branes may play a signihcant role in understanding this puzzle. 


2.4 Models encompassed 

In this section we will list the probe brane intersection models in which our results apply. 
We begin by recalling that in our notation (n | p T g) denotes the intersection of two stacks 
of Dp- and Dg-branes (with g > p) along n common directions. Let us also recall that the 
embeddings considered in this paper are the ones corresponding to massless quarks, in which 
the embedding function is trivial. We will hrst list the intersections which preserve some 
amount of supersymmetry. 


2.4.1 Supersymmetric intersections: 7)^ND=4 


The intersections (n | p T g) which preserve some amount of supersymmetry are those for 
which n is related to p and g as follows]^ 


n 


p + q — A 
2 


(2.44) 


^One can easily find this relation by imposing the no-force condition between the two stacks of D-branes 
(see, for example, 40 ). 
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It is straightforward to verify that the condition (2.44) selects the following three series of 
intersections: 


{p\p ± p + 4) , {p — l\p ± p + 2) , {p — 2\p ± p) . (2.45) 


Let us evaluate the index A for the three series of SUSY intersections (2.45). In these cases 
(1.1) simplihes drastically and we simply get: 


A = 2(g — n — 1) = q — p + 2 . (2.46) 

In other words A = 6,4,2 for the intersections Dp-D{p + 4), Dp-D(p + 2), and Dp-Dp, 
respectively, as illustrated in table 


Model 

A 

p 

q 

n 

Dp-D(p -I- 4) 

6 

p 

p + 4 

p 

Dp-D(p -I- 2) 

4 

p 

p + 2 

p — 1 

Dp-Dp 

2 

p 

p 

p — 2 


Table 1: The supersymmetric models for which our results apply. 


2.4.2 Non-supersymmetric intersections: 7 )^ND =6 


The models which break all the supersymmetries can be subject to instabilities. In flat space, 
the electromagnetic and gravitational forces of the two sets of Dp- and Dg-branes do not 
cancel out. This then typically manifests itself as a tachyonic mode below the Breitenlohner- 
Freedman bound in the open string spectrum in the near-horizon limit. However, certain 
circumstances may render the brane conhguration perturbatively stable; topology in the 
Sakai-Sugimoto model 41 or turning on an internal flux in the worldvolume of the probe 
The non-supersymmetric models where some of our results apply are listed 


D-branes [42,43 
in Table [2l 


Model 

A 

p 

q 

n 

Sakai-Sugimoto D4-D8/D8 41 

5 

4 

8 

3 

D3-D7’ 

42 


4 

3 

7 

2 

D2-D8’ 

44 


5 

2 

8 

2 


Table 2: The non-supersymmetric models for which some of our results apply. 

We do not plan to detail which of our results are directly applicable to the above non- 
supersymmetric models. This would invoke a separate involved study, so we just warn the 
reader by recalling a few facts. The Sakai-Sugimoto model is to be treated only in the 
deconhned parallel phase. The presence of the internal flux may affect the speed of zero 
sound and the attenuation; this expectation was recently conhrmed for the D3-D5 model 
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in 


45 . One should also keep in mind that all these models are subject to striped instability 


at non-zero density. Furthermore, in all the models, a further inclusion of the magnetic held 
will have a major effect due to Chern-Simons action contributions and hardly anything will 
apply. 


2.5 Bound on the speed of sound 


It was argued in 46 -48 that the speed of sound in a strongly-coupled theory with gravity 


dual is always smaller than the conformal value. In a (n -|- l)-dimensional theory this value 
is ul = 1/n. We will assume in what follows that n > 1. Let us explore under which 


circumstances our system violates this bound. Since we have obtained = 2/A (2.30), it is 
clear that the bound is violated if 

^<n . (2.47) 


From the general value of A written in (1.1), we conclude that (2.47) holds if: 


(p — 3)(p -I- g — 2n — 8) <0 


(2.48) 


In the SUSY case p + q — 2n — 9> = —4 (see eq. (2.44)). Therefore, the bound is always 


violated for the intersections listed in Table [T] with p > 3. For the non-SUSY intersections 
listed in Table only the Sakai-Sugimoto model violates the bound. The violation of the 
bound was also discussed in 37 for such intersections. 


Actually, one can probe in full generality that the bound is always violated for p > 
3 except for two particular intersections. In this case the condition (2.48) requires that 
p-|-g — 2U — 8<0. But, since the total number of spatial dimensions is 9, the integers p, g, 
and n must satisfy p-|-g — n — 9<0. Let us now write 


p-|-g — 2?7, — 8 = (p-|-g — u — 9)-l-(l — u) 


(2.49) 


The two numbers in parenthesis are less or equal to zero. Unless when they are simultane¬ 
ously zero, the bound is violated. This only happens when n = 1 and p -|- g = 10. The only 
intersections that satisfy these conditions for g > p > 3 are (1|4 T 6) and (1|5 T 5) for which 
cases ui = 1. 


3 Fluctuations 


We now want to analyze the excitation spectrum of our holographic system. These exci¬ 
tations correspond to density waves in the dual held theory which appear as poles of the 
retarded Green’s functions. In the holographic context hnding these poles is equivalent to 
obtaining the quasinormal modes of the gravitational system. Accordingly, we now assume 
that T and B are non-zero and allow huctuations of the gauge held along the Minkowski 
directions of the intersection, in the form: 

A = + a(p,x^) , (3.1) 
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where dx'^ = At dt + B dx'^ and a{p, x^) = ay{p, x^^)dx^. The total gauge held 

strength is: 

F = + / , (3.2) 


with F^^'> = dA'^^'> being the two-form written in (2.9). We will choose the gauge in which 


ttp = 0. Moreover, we will consider huctuation helds which depend on p, t, and x^. In 
this case it is possible to restrict to the case in which ^ 0 only when v = t,x^ = x, and 
x‘^ = y. It follows that the non-vanishing components of / are: 


ftp — Oi 5 fxp — Ctx 1 fyp ~ 1 

fxy dx Oiy , ftx 9 t CLx dx dt 5 fty dt dy , 


(3,3) 


where the prime denotes derivative with respect to p. In order to write down the Lagrangian 
for the huctuations, let us dehne the matrix X as: 


(p(°) + F(o)) V- 


Then, the DBI determinant can be expanded in powers of X as: 


y/-det(p + F) = A/-det(p(o)+ F(o)) 1 +-TrX - - TrX^ + - (TrX ) +0{X^) 
^ L 2 4 8 V / 

Let us split the inverse of the matrix -|- F^^^ as: 

-1 


(()<“>+ f<“>) =a-‘ + j, 


(3.4) 


(3.6) 


(3.6) 


where Q ^ is the symmetric part and JL is the antisymmetric part {Q is the so-called open 
string metric). It follows that: 

x\ = + T^fct, , (3.7) 

where the Latin indices take values in a, 6, c G {t, x, y, p}. The traces needed in the expansion 


(3.5) up to second order in X are: 

TrX = 


TrX^ = fed fab . 

In our case the relevant elements of the open string metric are: 


F~>tt ^ 

^ ” “7 

Jp 

gxX^ = g 


9r 


fp 9rr \9tt\ -^t 

9xx 


2 ’ 


g^^ = fp 


\9t 


9rr IS': 


tt\ 


A't^ ’ 


9L + ’ 

while those of the antisymmetric matrix ff are: 


j-tp _ _jpt — _ 


Ai 


9rr \9tt\ At 


12 ’ 


19 9 1 

J-X X _ _ ^X X 


B 


air + 


(3.8) 


(3,9) 


(3.10) 
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Let us write the components along t, p oi these matrices in terms of the charge density d. 
We get for the open string metric: 


= 


H + d? 


= ft 


H + df 


while takes the form 


j-tp _ _Jpi — 


9rrH ’ 

d VH + d^ 


(3.11) 


y/ 9rr \9 


tt\ 


H 


(3.12) 


From these values we can immediately calculate 


Ti-A' = 2 J: — kZLjF! + 2 , , 


\f9rr \9tt\ 


H 


9L + 


(3.13) 


and one can demonstrate that the linear term with TrW in the action does not contribute 
to the equations of motion, as it should. Moreover, it is straightforward to verify that, up 
to a multiplicative constant, the Lagrangian density for the fluctuations at second order in 
/ is given by: 


C 


\/9rr \9tt I 




rac rrbd 


red /raft 


The corresponding equation of motion for a'^ is: 
y/ 9rr I S': 


(3.14) 


dr 


Jtt\ 


fab 


= 0 . 


(3.15) 


From the equation of motion for Op (with Op = 0) we get the transversality condition: 


dta[ - u^{p) da;a'^ = 0 , 


where u{p) is the function 


u^ip) = 


9XX \9tt\ fp H 


g^^ gl, +52 H + d? 
Let us Fourier transform the gauge held to momentum space as: 

du dk 


ay{p,t,x) = 


( 2 vr) 


ay{p,oj, k) e 


—iu) t + ikx 


In momentum space the transversality condition (3.16) takes the form: 

oj a[ + vf{p)ka'^ = 0 . 

We now dehne the electric held E as the gauge-invariant combination: 

E = k at + UJ Qx ■ 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 
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(3.21) 


Using the transversality condition, we can obtain a[ and a'^ in terms of E' as follows: 
, kv? 


a* = 


E' 


/ ^ ^ 

" a;2 _ ^2 ^2 


E' . 


a;2 — /i)2 y2 

In terms of u the Lagrangian density of the flnctnations takes the form: 


C 


VSrr (glr + B^-)H 


u 


ac r*bd 




\/~ 9 xx~fp 

and the elements of the matrices Q and along tp are: 


jacjM ^ 


(3.22) 


= 


9lx + 


gpp = 


Qxx \9tt \ /, 


9rr{,9lx + B^) 


jj-ip gj^P^ — _ 9XX f 1 

“ \/9rr {9: 


f 

yxx 


+ B‘^)H 


(3.23) 


It is also interesting to write the explicit expression of the fnnction u for the Dp-brane 
backgronnd: 


u^ = 


9xx\9tt\fp H 

gl, + B^H + d? 


P fpip) 


p\ + pX+P-7 + d'^ • 


(3.24) 


The equations of motion derived from (3.22) for the Dp-brane background have been explic¬ 
itly written in appendix [A| (eqs. ( |A.4 ) and (A.6)). After hxing the Op = 0 gauge and using 
(3.21) these equations reduce to a system of two coupled second-order differential equations 
for E and ay which we study in the next three sections, both analytically and numerically. 

The numerical methods we employ to solving ordinary differential equations are by now 
standard. The only non-trivial complication comes from the fact that the fluctuations are 
generically coupled and one needs to find normalizable solutions for all the helds at once. Im¬ 
plementing such a method, though, is straightforward (see, e.g., 1^) and was hrst introduced 


in the holographic context in 49,50 


4 Vanishing magnetic field 


We will start our analysis of the fluctuation equations by considering the case in which 
the magnetic held vanishes. In this case the equations of motion for the longitudinal and 
transverse excitations decouple. Indeed, when B = 0 the momentum space equation (A.4) 
for the electric held E becomes: 


E" + dAog 


(p^ + d^)! /p 


(a;2 — fp kP) + uA 


E' + 


1 (a;2-/pfc2)p^ + n;2^2 




p^ + d2 


E = 0 . 


Similarly, when S = 0, the equation for ay written in (A.6) is: 
a" + dp log \\/p^ + d? 


a'y + 


1 (a;2 - fp k^) p^ + d^ 


P 




p^ + d2 


Op = 0 . 


(4.1) 


(4.2) 
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In the rest of this section we will analyze the solutions of (4.1) in different regimes. We will 


begin by analyzing in the next subsection the diffusive solutions of (4.1) in which u is purely 
imaginary and is related to the momentum k as u = —iDk"^, with D being the so-called 


diffusion constant. We leave for appendix the analysis of (4.2) and of the calculation of 
the corresponding transverse correlators. 


4.1 Diffusion constant 


Let us determine the diffusion constant which follows from the equation of motion of E at 


zero magnetic held (4.1). With this purpose we hrst expand this equation near the horizon 
p = Th- To begin with we expand the blackening factor fp{p) near p = r^. We get: 


/p(p) = (P - ^h) + 


Th 


(4.3) 


It follows that the coefficients of E' and E near p = can be represented as: 


dp log ■ 




^ — fp k?) p^ + d? 

1 (a;2 - fp p^ + (P 


P-Th 

A 


-|- Cl -)- ... 

C2 


/2 p^ + (P 


where the constants A, ci, and C 2 are given by: 

,2 


+ 


(p - ThY p-Th 


+ 


(4.4) 


A = 


CO 


(7-pYrl^ 

„A-1 


Cl = (7-p) 


k'^ p — 8 X r 


A-l 


r 


cP + 


p+A-6 


^ + 


+ w 


2rh 2 d? + rl 


C2 = 


h 


k^ + 


(7 _ p) (rf2 + ^A) - ■ {7-pYrl-^ 


u 


(4.5) 


Therefore, the near-horizon equation for E takes the form studied in appendix]^ (eq. ( |B.l )) 
and can be solved in a Frobenius series as in (B.2), he., as E{p) ~ (p — c/i)")! -|- (3{p — rP) + 


... ]. It follows from (B.4) that the solution of the indicial equation with infalling boundary 
condition is: 

a =-. (4.6) 


5-p 

2 


(7 - p) 

Let us next perform a low frequency expansion by considering k 0(e), a; ~ 0{e^). Then 
oi ~ C7(e^) , Cl ~ (!l(e ^) , C 2 ~ C^(c^) • (4-7) 

Since aci ~ C>(1), we get that: 




-a Cl . 


(4.8) 
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More explicitly: 




A+ 


p-7 


uj d? + 


(4.9) 


Notice that /3 ~ 0{1). Moreover, since a ~ a; ~ 0{e‘^), we can neglect the prefactor {p — rhY 
in the near-horizon expansion of E and write: 


E ^ Enh [l + Pip -rh)] , 
where Enh is the value of E at p = Vh- 


(4.10) 


Let us now perform the limits in the opposite order and expand the equation (4.1) for E 
in frequency hrst. One easily sees that the term without derivatives of E is of higher order 
in this expansion and can be neglected. Moreover, in the term with E' we approximate 
— fp /c^) p^ + u‘^ cP ^ fp p^ and then: 


log 


{p^ + d^)-2 fp 


-log 




(a;2 - fp p\ + ^2 5 

Therefore, the equation of E in this regime becomes: 


-1- constant 


(4.11) 


P 


E" - dp log 

{p^ -|- d^) 2 

For A > 2 this equation can be integrated as: 

poo 

Eip) = E^^^ + CE dp 


E' = 0 . 


(4.12) 


P^ 


3 5 

, 2 


(4.13) 


(p^ d'^) ■ 

where E^^'* = Eip = oo) (for A = 2 the previous integral is not convergent, see subsection 


4.2.2 for a detailed treatment for this). Next, we expand Eip) near the horizon: 

Eip) ^ + CEh 


Ce 


(r^ + d^) 


ip - rh) , 


(4.14) 


with 


h = 


P^ 


/ dp 


^ ‘-4^(3 l_i.3 _ l._^', 

i A-2 '* V2’2 A’2 ri / 


(4.15) 


By comparing (4.10) and (4.14) we get that — Eni, — Cp /a and that the integration 

constant Ce is given by: 


Ce = -P 



r 


A 

h 



P-7 

2 

h 



(4.16) 
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Thus, we can write the UV asymptotic value as: 


^2 


1 + i— (r^ + d^)^ h 

OJ 


By requiring that = oj^we hud the following dispersion relation 

OJ = —i D k"^ , 

where the diffusion constant D is given by: 


P-7 


D = r,- (r^+ /)»/,= 


2 P-5-A 

_ ly^ 2 

A-2 


/A ,0x1 /3 1 13 1 d'^\ 

(rh + d) ^(2’2 “ A’2 “ A’“^) 


(4.17) 


(4.18) 


(4.19) 


In order to write H in a more convenient way, let us dehne the rescaled density d as follows: 


d 


d = — = 
2 


7 — p\ 5^ d 


dvr 


A 

T ^~p 


Moreover, we also rescale the frequency and momentum as: 


OJ 


OJ = 


7 — p u 
“ An f ’ 


k = 


k 


7 — p k 


^ dvr T 


and dehne the rescaled diffusion constant D as the one which satishes: 

OJ = —i D k'^ . 

It follows that D and D are related as 

5^ dvr 

D = r.'^ D = -T D . 

^ 7-p 

Moreover, D only depends of d and A through the expression: 

2 ooxi /3 1 13 1 


D = 


A 


-2 ^ ^ V2’2 A’2 A’ J 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


We now analyze different limits of (4.24). First, we consider the case in which d —)■ 0, 


which is equivalent to the large temperature limit (see (4.20)). It follows immediately from 


(4.24) that: 


lim D = - - 

T-s-oo A — 2 


(4.25) 


Let us write this behavior in terms of the temperature T. Recall the relation between T and 


the ho rizon radius (2.18) (when the radius R is taken to be one). Therefore, it follows from 
(4.23) that the relation between D and D can be written in terms of the temperature as: 


7 — p '• 

D = -^D . 
AnT 


(4.26) 


^This is the widely used Dirichlet condition for the gauge field. We will consider other possibilities leading 
to anyonic correlations in section 
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Figure 1: The diffusion constant D in the D2-D4 model (A = 4, top) and D1-D5 model 
(A = 6, bottom) as a function of d. The red points represent the numerical data and the 
continuous curves are plotted using the formula (4.24). 


Thus, we have at large T that the diffusion constant behaves as: 

7-p 1 


D 


A - 2 27r T ’ 


T —)■ oo . 


(4.27) 


Let us now analyze the opposite regime in which d is large or, equivalently, when the tem¬ 
perature T is small for fixed charge density d. To obtain the behavior of the hypergeometric 


function in (4.24) in this regime, we use the following relation: 


F(a,l3-,r,z) = F(a,a+l-7ia + l-fth + 

1 (7^)1 (7 - a) ^ 

+ ~ a! HP, ^) + 1 - 7; /? + 1 - a; -) ■ (4.28) 

1 (a)i (7-p) 2 ; 


It follows that, for large d we approximately have: 


F 


/3 1 13 

1 

V2’2 A’2 

A’ 


3 1 


^ — T - -- T 1 + - d 


TT V2 A 


A 


z-i 


(4.29) 


and, therefore D takes the following approximate value for large d: 


^ (A-2)yF^(2“A)^(^^A)^' ^ (A-2) V2 A 


r(?-v')r(l + ^)d^r,-^ (4.30) 
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(4.31) 


P-5 


Since D = D and ~ T^-p, it follows that D behaves with the temperature as: 


_ 7 —p 

D T 


T ~ 0 . 


The numerical solution of (4.1) shows that the diffusion mode is the dominant one for high 


enough temper ature and allows the diffusion constant to be extracted. In £gure[T]we compare 


the prediction (4.24) to the results obtained numerically as a function on d for 


intersections. The numerical results for D agree very well with the prediction (4.24) which, 


wo different 


in particular, shows that the relation between D and d only depends on the index A. 

When the temperature is low enough, the dominant excitation mode is the zero sound. 
In the next subsection we will study analytically this mode following the method of . The 
transition between the collisionless and hydrodynamic regimes will be explored in subsec¬ 
tion |4]3l 


4.2 Zero sound 


The equation of motion of the electric held huctuation E at zero temperature and magnetic 
held can be obtained from (4.1) by taking fp = 1, namely: 


E" + dp\og 






+ d? 


(4.32) 


Let us study the solutions of (4.32) near the horizon p = 0 at low frequency and momentum 


(he., when a; ~ /c ~ 0{e)). Near the horizon p = 0, (4.32) can be approximated as: 

E” + a; ^ ~ 0 . 


(4.33) 


When p < 5, the two independent solutions of this equation are given by Hankel functions: 


E{p) = p^ 


(a) ( 2a; 

5^ V5 — p 


(a = 1, 2 , p < 5) . 


For p = 5 the two independent solutions of the diherential equation (4.33) are: 

= (P = 5). 


(4.34) 


(4.35) 


In what follows we will only consider the case p < 5. According to the standard prescription, 
the retarded Greens’ functions correspond to modes with incoming boundary conditions at 
the h orizon. Actually, in our case, when p —)■ 0 the argument of the Hankel functions in 
(4.34) grows (p^ —)■ oo for p < 5). As the asymptotic behavior of the Hankel functions is: 





^1— 
TIZ 


~ , (z^oo) , (4.36) 

V vrz 


it follows that the inco ming w ave corresponds to the function H^\ Thus, we select the 
solutions with a = 1 in (4.34). Moreover, when the index v of is not integer (in our 
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case this corresponds to z/ < 1 and p < 4), the Hankel function has the following expansion 
near the origin; 




1 

— + 


TT 


a\ 


2u 


71 a'' lx'' r(z/)r(z/ + l) V2/ 
for any constant a, when ax is small and z/ < 1. In our case: 


— ) [i — cot(7rz/)) x^' + 


(4.37) 


2uj 


a = 


P-5 

X = p ‘2 , 


u = 


5-p’ ^ ’ 5-p 

Therefore, after absorbing the multiplicative constant factor, we can write: 


E{p) = Ap + AcpU^-p + 


(p < 4) , 


(4.38) 


(4.39) 


where A is a constant and the coefficient Cp is: 


Cp — TT 


3-p 

(5 — p) 5-p 


ri ^ 

o—p 


i — cot 


71 


5 — p7 J 


(p < 4) 


(4.40) 


Notice that Cp has both real and imaginary parts. 


Let us next consider the case p = 4. The incoming solution contains in this case the 
function The corresponding expansion for small values of the argument of the function 

takes the form: 


H[^\ax) = —— - + ^ + 1 — 2'yE — 2 log {ax/2) + 


77a lx 


(4.41) 


where = .577 ■ ■ ■ is the Euler-Mascheroni constant. In our case x = p and a = 2u 
and we can write the expansion of E as 


E{p) = Ap + Ac 4 ^u‘^ + log + 


where C 4 is given by: 


C 4 = m + 1- 2'yE ■ 


(4.42) 

(4.43) 


Let us now perform the near-horizon and low frequency limits in the opposite order. For 
low frequencies and momentum, the equation for E becomes: 


E" + dp log 


(p^ d^) 2 


(a;2 — /c^) p^ -|- (P 
This equation can be readily integrated: 


E' TH Q . 


(4.44) 


E' 


(ca2 - k^) p^ + ^2 

CE -3- 

(p^ + d2)i 


(4.45) 
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where is a constant. A second integration gives: 


E{p) = - Cij 


(a;2 _ k^) pA + ^2 ^2 

{p^ + d?)^ 


dp , 


(4.46) 


where = E{p ^ oo) is the valne of the electric held at the UV bonndary. Let ns now 
dehne the following fnnctions: 


Mp) = 

Then it follows 


P 


Ip (p^ + d2)2 


dp , 


Mp) = 


dp 


Ip {p^ + d'^)^ 


E{p) = - CE - k‘^)Ji{p) + u^d^J2{p) 


For A > 2 these integrals are convergent and can be compnted analytically: 


(4.47) 


(4.48) 


Jiip) = 


p^-^E 


A-2 
2 


^ / X ^ 1 3 A „ / d 

1 hypergeometric fnnctions in 
{u^-e)J,{p) +a;2rf2^2(p) = -Y 


3 1 

1 3 

1 d^ \ 

2’ 2 

A’ 2 

A’ p^ J ’ 

/3 c 

1 5 

1 d? \ 

V2’ i 

A’ 2 

A’ p^) 


) J7i + cn 

^ d? J 2 can 

2 

P 

t-2 

A 

(p^ + d^) 2 


(4.49) 


1 A / 

' 1 1 

1 

3 

1 

A 

A-2 

. 2’ 2 

A’ 

2 

A’ 

p^) 


(4.50) 


Therefore, for A > 2, we have: 


E{p) = + 


2 ce 


P 


{p^ + d?) 2 




2k‘^ — Xu'^ 1 A / 

1 1 

1 

3 

1 

d^ \ 

A-2 

2’ 2 

A’ 

2 

A’ 

p^J 


(4.51) 


Let ns now expand E[p) in (4.51) near p 0. With this pnrpose it is better to deal directly 
with the dehnition (4.47) of the integrals J\ and J^- One can prove easily that: 

P 


Ji{p) = ^ ’ 


, . A — 2 ,l_q 
J2{p) = 


^ + 0{p^) , (4.52) 


where 7 is the constant dehned in (2.22). Plngging these eqnations into (4.48) we arrive at: 

,2 


E = CE^p + £("l + 


(4.53) 
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Moreover, taking into acconnt (2.21), we can rewrite E in terms of the chemical potential jj, 
as: 

. (4.54) 




Let us now match the two expressions we have found for E when p < 4 (eqs. (4.39) and 


(4.54)). From the terms linear in p we get the following relation between the constants A 

,2 


and ce'- 


ui 


A = ce^ . 

a 


(4.55) 


We can use this last relation to eliminate A. By comparing the constant terms in (4.39) and 
(|4.54l) we get 


2 ; 5 


^(0) = ^ 


d I X ^2 


2 ^~P 
+ CpU 5-P 


(p < 4, A > 2) 


(4.56) 


We now impose the Dirichlet boundary condition E^^^ = 0, which leads to the following 
relation: 

k-‘ - . (4.57) 

At lowest order, the right-hand-side of this equation can be neglected and we arrive at the 
following dispersion relation: 




(A > 2) . (4.58) 

Thus, the speed of zero sound coincides with the speed of hrst sound for the probe that 
we have determined in (2.30) and only depends on the index A of the intersection. Notice 
that the result (4.58) is also valid for p = 4 since the term on the right-hand-side of (4.57) 
is higher order. To avoid clutter, let us consider the solution with plus sign in (4.58) and 
expand u as: 

oj = \j~^^ + ■ (4.59) 

Plugging this ansatz in ( |4.57[ ) and expanding at zero order in 5a;, we arrive at an equation 
for 6oj\ 

V^A Cp 

~ ' '.A. 


du = 


2\f2 p 


\ 5—p 1 '^~P U p / 2i \ 2(5—p) , ^ 


fcs-p -|- 2 


5 — p VA 


k^-p 6u 


(4.60) 


which we can solve for 5a; in powers of k. The second term on the right-hand-side is clearly 
of higher order, thus 


6u = —^ ( Y 

2/i VA 


O \ '^~P _ 

2(5-p) , IIZP 

k^-p 


(4.61) 


Taking into account the expression of Cp in (4.40), we hnd the following expression for the 
imaginary part of uj: 


3-p 


Ima. = 0 + Imfo = -,7A^, 

2p -p ( 1 ^ ^ 

' 5-p 


O \ '^~P „ 

^ \ 2{5-p) , ^ 
k^-p 


(4.62) 
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Similarly, we get a correction to the real part of w, namely: 


3-p 

TT (5—/ TT 
Reoa; = -- ft cot 


2p 


r( ^ 

0—p 


5 — pJ VA 


7-p ^ 

2(5-p) , l=P 
fcs-p 


(4.63) 


In Fig. [^we plot our prediction for the real part of u and we compare it with the results 


Re w Re cj 




Figure 2: In this plot we represent the real part of Cj versus k for two classes of systems. On 
the left we depict the dispersion relation for the intersections D3-D3 (p = 3, A = 2), D3-D5 
(p = 3, A = 4), and D2-D6 (p = 2, A = 6). On the right we plot the same quantity for 
D4-D4 (p = 4, A = 2), D4-D6 (p = 4, A = 4), and D4-D8 (p = 4, A = 6). The models are 
listed top-down. The points correspond to numerical values at d = and the continuous 

curves are the analytical results. 

found numerically for different values of k. We notice that the numerical dispersion relation 
for the zero sound mode is well reproduced by our analytic equations. 


4.2.1 The p = 4 case 


When p = 4 the general formula (4.54) does not match (4.42), due to the presence of the 


by: 


logarithmic term in the latter (which is absent in (4.54)). In order to include these terms we 
have to compute the next order correction to (4.54) near the horizon, following the procedure 
of 


12 . First of all, we recall that the near-horizon equation (4.33) for E for p = 4 is given 


ui 


E" = - -E . 

p3 


(4.64) 


Neglecting the right-hand side in this equation and integrating it gives rise to a linear solution 
as in (4.54), which is second order in the expansion parameter e. To go beyond this order 


we just take into account the right-hand side of (4.64) and substitute the value of E, taken 


from (4.54), in this term. We arrive at the equation: 


E" — —up' [Aq p ^ Ai p 


(4.65) 
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where Aq and Ai are: 


A„ = £!») + 

X a 2 ■ 


u 


Ai — Ce —r 

d 


(4.66) 


(Aq + Ai p is just the second-order expression (4.54)). Integrating this expression twice, we 
get 

m2 

(4.67) 


^ . 9 . Aq 

E = Ao + Aip -h Ai (n log p---h ... , 


2 p 


where we took into accou nt tha t the solution of the homogeneous equation is just (4.54). Re¬ 
call also that the solution (4.64) with in-falling boundary conditions is given by p^ /p^) 

and that the expansion (4.42) is valid when the argument of the Hankel function is small, i. e., 
when p 3> Therefore, we can neglect the last term in (4.67). The resulting expression 
can be written as: 




2ce P 


A 


OJ 


- ^^"]+CE^logp + 


(4.68) 


Let us now match (4.68) and (4.42). By comparing the terms linear in p and those logarithmic 
in p, we hnd that the constant A is related to cg as in (4.55). Moreover, from the identihcation 
of the constant terms, we get that is given by: 


Eio) ^ ^ 

d 


2p ,X 


A '2^ 


+ (c4 — logcn^) 0;“^ 


(p = 4,A>2), (4.69) 


where C4 has been written in (4.43). By imposing that E^^'^ = 0 we get the following 
dispersion relation: 


= — (c4 - loga;^)^'^ , (p = 4,A>2). (4.70) 

Let us solve this polynomial equation at different orders in the expansion parameter e. At 
leading order we hnd that uj is real and given by (4.58), he., this general expression is also 

(4.71) 


valid for p = 4. The next-to-leading contribution 5u is given by: 


(5a; =-- 

2p vA 


2\f 


C4 


>og(^))^=. 


(p = 4, A > 2) . 


Separating the real and imaginary parts using (4.43), we get: 


TT / 2 

Im a; = 0 -I- Im (5a; =-— 

2p \X 




Re (5a; = 


27 e - l-hlog(2P/A) /2\f 


2p 




(p = 4,A>2) 


(4.72) 


Notice that, curiously, the leading term to Imo; written above is exactly the same as the one 
obtained by taking p = 4 in (4.62). For A = 5 these formulas were obtained in (l^. 
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4.2.2 The A = 2 case 


When A = 2 the integral J7i(p) defined in (4.47) is not convergent and, therefore, the integral 
of (4.45) written in (4.48) is not valid. In this case we can define a new integral Si{p) as: 


Mp) = 






(p2 + d2)i p\ ^p2 + ^2 


P 


1 + log 


2p 


\/p^ + (p + p 


(4.73) 


Then, the solution of (4.45) for A = 2 can be written as: 


E{p) = - ce - k‘^) (Ji(p) - logp) + (f J 2 {p) 


(4.74) 


At the UV (he., when p is large) the integrals J7i(p) and 772 (p) vanish by construction 
and E{p) therefore behaves as: 


E{p) = E^°'^ + Ce (w^ - fc^) logp + 


(p ^ oo) 


(4.75) 


When this logarithmic behavior is present, the source in the AdS/CFT correspondence is 
determined by the coefficient of the logarithm, which should vanish. This can be achieved 
either by requiring Ce = 0 or: 

u = ±k . (4.76) 

Let us see that if for ce = 0 the solution of the fluctuation equation is trivial. Indeed, in 


this case E{p) is constant and the matching with the near-horizon results (4.39) or (4.42) 
is only possible if the constant A in those solutions is zero. This, in turn, implies that the 
whole solution vanishes, as claimed. Therefore, the dispersion relation in this case should be 


given by (4.76), which corresponds to having no dissipation and the speed of zero sound is 


equal to one. Notice that, when (4.76) holds, E{p) behaves near p = 0 as: 


E{p) = +ce‘^P + 0{p^) . 
a 


(4.77) 


One can then match (4.77) with eq. (4.39) (for p < 4). Fr om t he comparison of the linear 
terms we get that A is given by the same expression as in (4.55) and that E^^'^ is related to 
Ce as: 


d 


2 ^~P 
OJ 5-P 


(4.78) 


Notice that (4.78) is just the same as (4.56) when A = 2 and oj = ±k. However, in the 


present A = 2 case we do not have to impose that E^^'> = 0 and, thus, the dispersion relation 


is not given by (4.57) but instead we have (4.76) 
intersection). 


see 


13 for a similar analysis in the D3-D3 


When p = 4 the matching of the logarithmic term in the near-horizon expansion requires 


to go beyond the leading term in oj, as in subsection 4.2.1 It is easy to see that (4.77) can 


be corrected in such a way that it matches (4.42) and that E^^'^ is related to ce as in (4.69) 
with A = 2 and u = Ek 
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The dispersion relation (4.76) is the natural one for a massless excitation in 1 + 1 dimen¬ 
sions. Accordingly, let us determine all possible intersections with A = 2 and n = 1. Taking 
A = 2 and n = 1 in (1.1) we hnd the following relation between p and q: 


{p-3)q = (p-3)(10-p) , (A = 2,n = l) . (4.79) 

When p = 3 this equation is satished by any g, whereas for p 7 ^ 3 it requires that p + q = 10. 
Therefore, we have the following two series of intersections having A = 2 and n = 1: 


(l|3Tg), g = 3,5,7. 


(l|pT(10-p)) , 


where we have required that q > p. 


P=l,2,4, 


(4.80) 


4.3 The collisionless/hydrodynamic crossover 



k 


- 0.008 


Figure 3: We are presenting typical dispersions in the D2-D6 model at d = 10® and 5 = 0. 
The dashed curves are real parts and continuous curves are imaginary parts. Around k ~ 
0.006 the diffusive mode merges with another purely imaginary mode, which dehnes what is 
called the crossover location {kcr,i^cr)- Notice, that we have included the first excited purely 
imaginary mode from the transverse sector (lowest pink curve). 


Our analytic treatment of the zero sound in subsection 4.2 was carried out at T = 0. The 
numerical analysis of the fluctuation equation shows that the zero sound persists at T 7 ^ 0 
if T is low enough. A typical dispersion relation is depicted in Fig. In this regime the 
real part of this mode is independent of the temperature, while its imaginary part receives 

7-p PI 

corrections which are proportional to T^-p. This behavior is illustrated in Fig. for the 
D1-D5 intersection. 


If we continue increasing the temperature, at some point there is going to be a crossover 
to a hydrodynamic diffusive regime, in which the D behaves as in eqs. (4.31) and (4.27). 













Im[ 6(j] 


Log[ -] 

lm[ w ( 7 = 0)] 



-- 2/ A 
Log [ d ] 


Figure 4: The imaginary part of the lowest excitation mode in the D1-D5 model at fixed 
= 0.01. The red point represent the numerical results and the line segments have 


(i(5-p)/A 

slopes 0,3,-3,-2 from left to right. These are consistent with the analytic results. The 0 
conforms with the notion that the zero sound is robust against temperature variations and 
only for sufficiently hot environment will the decay rate behave as The latter two slopes 
correspond to the diffusive regime with temperature behaviors D T (eq. (|4.3l|)) and 
T-i (eq. (|T27|)). 


D 


The frequency Ucr and momentum kcr at which this collisionless/hydrodynamic crossover 
takes place depends on the temperature and chemical potential. We have determined this 
dependence numerically (see Fig. |^. From this analysis we conclude that ujcr and kcr scale 


as: 


7-P 
T 5-P 


Uc 


kcr ~ 


7-P 

T 5-p 


(4.81) 

r p 

Notice that, for the D3-D5 and D3-D7 system, Ucr ~ in agreement with the analysis 

respectively, and for the Sakai-Sugimoto model Ucr rVp. in agreement 


of 19 and 


22 


with the analysis of 29 . 


5 Non-vanishing B field 


In this section we study the influence of the magnetic field B in the collective excitations 
of our brane intersection. The equations of motion of the brane fluctuations have been 
written in appendix [A| (eqs. (A.4) and (A. 6 )). Notice that, when the magnetic held is 
non-vanishing, the longitudinal and transverse huctuations are coupled. We hrst study the 
magnetized system at T 7 ^ 0 in the diffusive regime. 


29 
















p 



^-2M 

Log d 


- D1-D5 - D2-D4 - D3-D5 - D4-D6 
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Figure 5: We plot the frequency (top panel) and the momentum (bottom panel) for the 
location of the transition from the hydrodynamic to collisionless regime. The scalings are 
chosen such that the slopes for all the models is 1 (best £t values). The models are BI¬ 
DS, D2-D6, D3-D5, and D4-D6 top-down. This translates to scalings uicr ~ ~ 

H and kcr ~ ~ /jj,. Notice that we have subtracted a 

constant p to separate the lines. 

5.1 Diffusion constant 


To obtain the diffusion constant in the presence of the magnetic held we follow the same 
steps as in subsection 4T First we expand the equations of motion around the horizon 
p = Th. As in ref. 22 the equations decouple and we have to expand the terms multiplying 
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E' and E in (A.4). These expansions are: 

„7-p 


dp\og 


P 


(pA^^A+p-7^2 
p7-p + 52 (^2 _ ^2) pX + ^2 52 pA+p-7 + ^2 ^2 


1 (a ;2 - fp k^) + a ;2 52 ^A+p-7 ^ ^2 p2 


A 


p’^~^ fp P^ + -02 + ^2 

where the constants A, ci, and C 2 are given by: 


+ 


P-Th 

C 2 


+ Cl + 


[p-rhY p-Th 


+ 


(5.1) 


A = 


UJ 


{7-pfrl-^ 
Cl = (7-p) 


„A-1 


A :2 

-77 + X- 


p — 8 + (6 — p)B‘^rf^ 


d2 + + 52 a;2 2r;, (1 + 52 { d ^ + 

p — 8 + A + (A — 1) 52 


d 2 + + 52 


„P+A -6 


C 2 = 


^ ' h 

.2 


(7 - p) (d2 + + 52 


'fc 2 + 


(7-p)2rr^ 


o; 


( 6 . 2 ) 


The resnlting flnctnation eqnations can be solved in Frobenius series as in (B.2), in terms of 


two parameters a and (3. In the low freqnency limit with k 0{e),u: C>(e2), the exponent 


a is given by the same expression as in (4.6), while (3 is given by: 

A +2 


h 2 

(3 = 1 — 


p-7 


00 + r^ + 52rr^-^ 


(5.3) 


We now expand hrst in freqnency. The eqnation for E deconples from the one for ay in 
this limit and becomes: 


E" - dp log 
This eqnation can be integrated as: 


E' = Q . 


E{p) = + CE 


pA + pP+A-7 52 

f {p^ + pP+^-'^ 52 _|_ g^ 2 ^ I. 

pA ^+A-7 j^2 


dp 


[p^+pp+^-7 B^+d^y 


(5.4) 


(5.5) 


where ce is a constant of integration. Let ns expand the fnnction E{p) in (5.5) near the 
horizon as: 


^A I A+p-7 o2 

E{p) = + ceIx,p - Ce , , (p - rh) + 


(r^ + r^+^-"52 +d2)f 


(5.6) 


where Z\^p is the integral: 

pA + pP+^-7 52 


Xa,p = / dp — -^^ 

Jrh (^pA -)- pP+A-7 J^2 _j_ d"^^ 2 


dx 


X 


A ^ ^p+A-7 q2 


11 _|_ 2,p+A-7 ^2 _|_ ^2 j 


3 5 
2 


(5.7) 
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and where we have rescaled the magnetic held B as: 


B = 


B 
. 2 


7-p\^ B 
Att 


, ,7-p 


(5.8) 


We now match (5.6) with the near-horizon expression (4.10). By comparing the constant 
terms in both expressions we conclude that = E^h — ceIx,p, whereas the constant ce 
can be determined by looking at the linear term. We get: 


Ce = -13 E, 


nh 


(1 -I- -I- d^)2 

1 + 52 


(5.9) 


where f3 is given in (5.3). It follows that: 


= K 


nh 


J .2 (1 + 52 ^ 


1 -|- i — r. ^ 
u 


1 + 52 


-\,p 


(5.10) 


Requiring E^^'> = 0 we h nd a dispersion relation of diffusive type, a) = —iDx^p K, where a) 
and k are rescaled as in (4.21) and the rescaled diffusion constant Dx^p is given by: 


D\,p — 


{1 + B^ + d^y 

1 + 52 


dx 


X 


A ^ ^p+A-7 q 2 


11 ^^A _j_ ^p+A-7 ^2 


(5.11) 


We have not been able to compute the integral ( |5.11 ) analytically for arbitrary values of A 
and p. However, when A -|- p = 7 this integral can be obtained in terms of a hypergeometric 
function: 


D\=7—p,p 


2 {l + B^ + iy 


1 + R 2 
2 

+ 




(d2 + 52)(l + d2 + R2)2 


R 2 


A-2 d^ + B"^ V2’2 A’2 A 


11 13 


-d^ - B^ 


(5.12) 


Notice that the systems for which (5.12) holds include the D3-D3 intersection ((2|3 _L 5)) 
and D2-D8’ ((2|2 ± 8)). In Fig. |^we compare the analytic results to the numerical values 
of the diffusion constant as a function of B. It is also easy to hnd the limiting value of Dx^p 
as d = 0. Indeed, for any value of A and p we get: 


Dx,p{d = 0 ) = 


A-2 


A 1 + 


A-2 


-^— e(- 

V2’2(7-p)’ 2(7-p)’ 


-B^ 


(5.13) 


Notice that when R = 0 eq. (5.13) reduces to (4.25). 

We wish to end this subsection with the following remark. In analogy to the lower bound 
on the shear viscocity to the entropy ratio, Hartnoll proposed a universal lower bound on 
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Figure 6: Plot of the diffusion constant for D2-D4 and D2-D6 as a function of B at fixed 
d = 10. The continuous curves are the analytic results, whereas the points are from the 
numerics. 


the charge diffusivity for incoherent metals at high temperatures 51 . Unfortunately, the 


proposal was not precisely formulated and thus the bound does not have a specihc value, and 
we therefore cannot make a direct comparison. Nevertheless, according to 
in our language reads 


51 , the bound 


D > 


hvp 

kn 


= vl 


(5.14) 


where vp is some number of order 1 (for Fermi liquid with quasiparticle description vp would 
correspond to the Fermi velocity). The minimum value for D occu rs at d = 0, i.e., at high 
temperature as in the regime of interest. This was computed in (4.25): lim-r^oo d) = 2/(A—2). 
At hnite magnetic held this is given above in (5.13). In general, turning on the magnetic 
held the dihusion constant decreases without bound, see Fig. for example. Therefore, 
holographic metals as modelled in the present work, seem to evade any bound on charge 
dihusion at least at hnite magnetic held. 


5.2 Zero sound with B field 

Let us study the zero sound for intersections with B held. Typical dispersions are depicted 
in Fig. [^for the D2-D6 model. First of all we analyze the equations of motion of appendix 
[^for E and ay at zero temperature near the horizon p ^ 0. We will assume that B is small 
(p7-p ^ these assumptions the equations of motion for E and ay greatly simplify 
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Figure 7: The zero sound becomes massive above a critical magnetic field strength. We 
depict the dispersions for the D2-D6 model aX B = 0.1 (left panel) and B = 0.2 (right 
panel), both a.t d = 10®. In the latter case the zero sound mode has just become massive 
and starts propagating. 


and become the coupled system: 

{7-p)B^ 


E” + 


p{p'^~P + B"^) 


E' + 


u 


P 

{7-p)B‘^ , 


7-p 


E = —i{7 — p) 


Bu^ 

p(p7-p + 52) 

B 


ci.j ~cb., + — (ly — i{7 — p) —. 

+ B^) y p^-p ^ ^ ' p{p^~y + B^) 


E . 


Let us now define the operator O, which acts on any function E[p) as follows: 

dE{p) = E" + 


{7-p)B- u 


p{p'^~P + B'^) 


2 , ,2 
E' + ^ p 


P 


7-p 


Then, the equations in (5.15) can be written as: 

Boj^ 


OE = -i{7-p) 


p{p'^~P + B‘^) 




O tty = i{7 — p) 


B 


p{p'^~P + B"^) 


(5.15) 


(5.16) 


E . 


(5.17) 


These coupled equations are easy to decouple. By defining y± as 


E 

2 /± = — ± Op , 

lOJ 


the equations satisfied by y± are decoupled, and given by: 


C>± (7-p) 


Boj 


p{p'^~P + B'^) 


y± = 0 


(5.18) 


(5.19) 


More explicitly, we have 


y± + 


{7 -p)B^ 

p(p7-P + 52) 


y'± + 


^ (7 — p) Boj 
p'^-P pip'^-P + 52 ) 


h/± = 0 . 


(5.20) 
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We now study the equations for y± when the B held is small. Let us neglect the terms that 
are quadratic in B and approximate these equations as: 


!/^ + (5.21) 

Let us expand the solutions of these equations in powers of B. Notice that the equation for 
H- is obtained from the one for by changing B by -B. Accordingly, we write: 

y±{p) = yo{p)±Byi{p) . (5.22) 

The equations for yo and yi are independent of B and given by: 

2 2 

y'o + ^2/0 = 0 , y'i + ^ 1/1 = -(7-p)^|/o . (5.23) 


The solution for the equation of j/o with incoming boundary condition is just the Hankel 


function written in (4.34): 


- tj{i) 
yo = c+p^ 

5-p 


2u 

5 — p 


P-5 

P " 


(5.24) 


where c+ is a constant. Notice that yo is a source in the equation of yi. Actually, dehning 
the function y as j/i = c+j/, we hnd the following inhomogeneous equation for y{p)'. 


y 


/_ N 1 tt(i) ( p-5 \ 


(5.25) 


In the appendix we hnd the solution of this equation by using the Wronskian method. 
This solution is simply: 


y = 


£21® rr{l) 

-p 2 H\^ 

5-p 




(5.26) 


which also satishes the incoming boundary condition at the horizon. Let us now write the 
solution for y±. First we dehne the functions: 


Zi{p) - p^H^} ( P ^ , 

5^ \0 — p / 

z,(p) = p’-^) . 

522^ V5 — P / 

(5.27) 

Then, y±(p) are given by: 



y+{p) = c+zi{p) - c+Bz 2 {p) , 

y-{p) = c_zi{p) -h c_Bz 2 {p) , 

(5.28) 


where, to obtain y_ we changed B ——B and c_|_ —)■ c_. Let us now redehne these constants 
as follows 

ibJ ibJ 

Cl = y (c+ + C_) , C 2 = (c+ - c_) . (5.29) 

Then, E and ay can be written in matrix form as: 



^i(p) 

W^ 2 (P) 


^2^2(p)\ 

L-i(d)/ 



(5.30) 
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J, 

We now study the near-horizon solution found in the small a; limit. At 0{u^) we have 
(for p < 4): 


, ^ 2 "r(z/) . 

a(p) ~-* 


TT 0 !“ 


[p + CpWS-p] , 


Z2{p) 


2^T{u) . 1 

■ % 


TT a'' 


OJ 


(6.31) 


where u and a are dehned in (4.38) and Cp is given in (4.40). Then, after redehning the 


constants C\ and C 2 by absorbing the common factor in zi and Z 2 , we have: 

' E\ (p + Cr,uj^ ^ ^ 


V 


iB 

„i2 


h {p + cpu;^-p)j 



(5.32) 


We will now perform the two limiting operations in the opposite order. As in 22 , in the 
low-frequency limit we drop all terms not containing derivatives of E and ay, as well as all 


factors of B. Thus, the equation of motion of E reduces to (4.44), whose solution is (4.51). 
The corresponding equation for ay in this limit becomes: 


a'' + dp log [p^ + d 


ay = af - 


This equation can be integrated twice to give: 

Cj/ jTs (p) ) 

where Cy is a constant and Jz{p) is the following integral: 

dp 



1 

3 

1 

d^ \ 

A- 2 ^ ^l 2’2 

A’ 

2 

A’ 

p^J 


It is clear from these equations that a^^ is the value of at p —)■ cx). Moreover: 


(5.33) 


(5.34) 


(5.35) 


ay{p) = a. 


_ ( 0 ) _ 2 Cp 


/I 1 1 3 


A-2^ '^V2’2 A’2 


1 


(5.36) 


Let us now expand ay{p) in powers of p. First, one can check that, for small p, the integral 
5 / 3 (p) can be approximated as: 

(5.37) 




where p is the chemical potential ( |2.21| ). Therefore, for small p: 

ay{p) ~ 4 °^ 


p , p 
""" d +""" d • 


(5.38) 


The corresponding expansion of E has been written in (4.54). We now match (4.54) and 


(5.38) with (5.32). By comparing the terms linear in p we hnd the constants c\ and C 2 in 
terms of Ce and Cy, namely: 


u 


Cl , Ce , 

d 


C2 ^ 7 Cy , 


UJ 

d 


(5.39) 
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and the comparison of the terms independent of the radial variable leads to the following 
matrix equation: 




.( 0 ) 


/ 


2(6-p) 






ill 

d 


\ 


+ -dj 


Ce 


(5.40) 


We now require the vanishing of the sources and a^\ which only happens non-trivially 
if the matrix written above has a vanishing determinant or, equivalently, when: 


2(7-p) 2(6 -p) 

CpU 5-P + 2ldCpU 5-P — 


^5-p _|_ ^2^2 _ ^ ^2 ^2 _ 2^2 ^ Q _ (5.41) 

A A 


Eq. (5.41) is a polynomial relation in u and k whose solutions determine the dispersion 


relation of the zero-sound. Let us solve this equation for u. At leading order we can neglect 
all the terms except the last three, which leads to the following gapped dispersion relation: 




A 




(5.42) 


In order to solve (5.41) at the next-to-leading order, let us write: 


/ 2 ,o ^ 

u — \ ~r k"^ -|- — r- -|- ou! . 


A 




(5.43) 


By plugging (5.43) into (5.41), we get at next-to-leading order 



B2^ 

2(5-p) 


Bh 

,2 

\ Tk + 



— 

—— 



/ 

- A 

p2j 


(5.44) 


Let us now separate the real and imaginary parts in Su. From the expression of the constant 
Cp in (4.40) we obtain that the imaginary part at this order is given by: 


3-p 


Imca =- 

h 


r( ^ 

0—p 


h .2 

B2^ 

2(5-p) 

ik^ 

Bh 

\ Tk + 



— 

- — 



/ 

- A 

p2j 


Moreover, we get the following correction to Reca: 


3-p 


TT (5 — p)5-P / TT 

Redo; =-^ cot 




r( ^ 

5—p 


5 — pJ \ X 






L A ^ /i2 


(5.45) 


(5.46) 


Notice that (5.45) and (5.46) coincide with (4.62) and (4.63) when the B vanishes, as it 
should. 
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Figure 8: We present the real part of the dispersions in the D2-D6 above the critical value 
above which the zero sound mode has become massive {13 ~ 0.2) at various B = 0.5,5,10 
(bottom-up) for d = 10®. The points are the numerical values and the continuous curves are 
the analytic results. Notice that the lowest curve has a small dip at small k diverging from 
the analytic result. This is easily understood by recalling that the zero sound picked up the 
mass just slightly earlier. 

5.3 The crossover with B field 


Now let us hnish this section with the following observation. Recall that in section 4.3 


we were investigating the location where the transition from the collisionless regime to the 
hydrodynamic regime takes place. By studying the dispersion relations in the complex c5 
plane by increasing the real momentum k we found that the diffusion mode met with another 
purely imaginary mode at some non-zero ((her, kcr), where they merged together to become 
the pair of zero sound modes propagating in opposite directions. 

While we now know from several studies that the zero sound generically picks up a mass 
with a non-zero magnetic held, it is interesting to ask if the location of the transition point 
is affected by the magnetic held, too. It turns out that the answer to this question, properly 
dehned, is negative. By increasing the magnetic above a critical one, where the zero sound 
mode has become massive, the dihusive mode never meets with another purely imaginary 
mode (see Figs.[^&:[^) and it is unclear how to dehne the location for the crossover to happen. 


However, as suggested in 


22 


one can dehne it by making the momentum k complex and 


keeping Cj real. Indeed, we repeated a similar analysis to theirs (D3-D5 and D3-D7 cases) 
in the D2-D6 model and found out that the transition point is unahected by the magnetic 
held strength. This suggests that this phenomenon is rather generic and may happen for all 
other intersections as well. We have depicted the dispersions in the D2-D6 case in Fig. 


38 








Log lm[/(] 



Figure 9: We are depicting the dispersions as complex momentum as a function of frequency 
for the D2-D6 model. On the left panel we plot the real part of k and on the right Im/c. 
The red, lower, points are numerical data for B = 0 and the blue, upper, points are for 
B = 10. Both cases are for d = 10^. The continuous lines starting from the left are analytic 


results for the diffusion mode (eq. (5.11)) and the continuous curves starting from the right 
are the analytic results (at T = 0) for the zero sound (eq. (5.43)). The lower point of 


the vertical line segment is extracted from the collision of poles in the complex u plane, 
indicating the collisionless/hydrodynamic crossover location. We infer that the transition 
point is essentially insensitive to the inclusion of the magnetic held. 


6 Alternative quantization 


In this section we will not allow the intersection of the probe with the background branes be 
as general as before, but we will specialize to a (2 + l)-dimensional intersection. Contrary 
to the scalar helds, for the gauge helds one can impose mixed Dirichlet-Neumann boundary 
conditions, i.e., adopt an alternative quantization, only when the bulk spacetime is four¬ 
dimensional 52,53 . From the boundary perspective, a mixed boundary condition for the 


gauge held can be traced back to a SL{2,Z) mapping from the usual Dirichlet condition, 
where this operation generically corresponds to adding a Chern-Simons term to the action 
and making the external vector held dynamical]^ We refrain from presenting a lengthy 
introduction to the details (see also |34| - [3^ ). Instead, we will walk the reader through 
the necessary notations along way as our goal is a simple generalization of the methods 


developed in 33 to tackle all conformally AdS 4 ^ backgrounds. 


We wish to emphasize the key observation made in (allowing the use of alternative 
quantization) that the bulk action itself does not have to be invariant under SL{2,Z) trans¬ 
formation: one only needs the bulk equations of motion for the gauge helds to reduce to free 
equations near the boundary so that one can choose mixed boundary conditions. In other 
words, since diherent quantizations only diher by boundary terms, the equations of motion 


(A.4) and (A.6 ) which we plan to solve are still the same. The diherence comes in interpret¬ 


ing the diherent quantities: in general the quantum liquid transforms into an anyonic one at 
non-zero density and magnetic held. For precise mapping between the parameters, we refer 
the reader to 


33 


^Incorporating SL{2,Z) duality has sparked many recent studies in holography [33|j3^[54||^ . 
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Let us thus consider fluctuation modes satisfying the following mixed Dirichlet-Neumann 
boundary conditions: 

( 6 . 1 ) 


lim 

p^OO 


/pM - X = 0 1 


where n is some constant. The Dirichlet boundary conditions considered so far in the normal 
quantization correspond to n = 0. Equivalently, as in 35 , we can write the previous condition 
as: 


lim [ fp^ - ^ ] = 0 . 

p^oo L Z 


( 6 . 2 ) 


Clearly m = 1/n and the normal quantization corresponds to m = cx). For /j, = t,x the 
condition (6.1[) leads to: 


lim \n a[ — i k tty] =0, 

p^OO ^ 


while for p = y it becomes: 


Let us rewrite (6.4) as: 


lim \n P '2 a'^ + iuj Qy] = 0 , 

p^oo 

lim \np^ a — iE] = 0 . 
lim E = —in lim [] • 


p^OO 


p^oo 


(6.3) 

(6.4) 

(6.5) 


Moreover, by using the transversality condition (3.21) at p —>■ oo, where u —)■ 1, the two 
equations in (6.3) reduce to: 


lim tty = i 

p^OO 


OJ 


^ p-l-oo 


lim [p2 E'] . 


( 6 . 6 ) 


As stated above, it is clear from (6.5) and (6.6) that n = 0 corresponds to the Dirichlet 
boundary condition in which E and ay are fixed at the boundary. 


Recall that the horizon radius can be eliminated from the fluctuation equations (4.1) 
and (4.2) by rescaling the radial variable p as p = p/rh and rescaling d, oj and k as in (4.20) 
and (4.21). We will rescale the gauge potentials as dp = ap/rt, which for the electric held E 
means the following: 

E 


E = 


7-p 
. 2 


(6.7) 


It can be easily checked that we can eliminate from the boundary conditions (6.5) and 

n 


(6.6) by rescaling n as: 


n = 


7-p-\ 
. 2 


h 


Equivalently, m must be rescaled as: 


m 


m = 


P+A-7 
, 2 


( 6 . 8 ) 


(6.9) 
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Let us now analyze the new boundary conditions at low frequency and momentum. The 
expressions of E and ay in this regime have been written in (4.46) and (5.34). The UV 
values at the boundary of these two helds are: 


lim E{p) = , 

p^OO 


lim ay{p) = a® 


p^OO 


( 6 . 10 ) 


The radial derivatives of E and ay at the boundary can be easily obtained from (4.45) and 

- k^)cEp~^ , 


(p4|: 


dE 

dp 


da. 


Cyp 2 


p^OO 


( 6 . 11 ) 


p—>-00 op 

From these expressions we can recast the boundary conditions for the alternative quantization 
as a relation between the constants E‘'^\ a^\ ce, and Cy. Indeed, let us dehne eP and a^yJ^ 
as: 

incE ■ ( 6 - 12 ) 


eP =E^^^ + tr\cy , = a® 


Then, (6.5) and ( 6 . 6 ) are equivalent to the conditions: 

— n 


(6.13) 


6.1 Zero sound 


From the results of section 5.2 it is straightforward to relate E^ and a'^yj^ to the constants 
ce and Cy for incoming boundary conditions at the horizon at low u and k. Indeed, from 


.( 0 ) 


(5.40) we get: 




/d“’\ 

2 (6 -p) 



\aS/ 

V 

f-in 

+ -J 


Ce 


(6.14) 


The alternative quantization conditions (6.13) can be implemented by imposing the vanishing 


of the determinant of the matrix on the right-hand side of (6.14). The corresponding equation 


is just (5.41) with the substitution B ^ B — nd. Solving this equation we get the dispersion 


relation. At leading order the real part of u is just: 


UJ = ±-1 


, (B-ndr 




(6.15) 


Similarly, we could get the higher order terms from (5.45) and (5.46). Notice that the 


spectrum can be made gapless by adjusting the alternative quantization parameter to cancel 


the gap induced by the magnetic held. See Fig. 10 for the comparison between (6.15) with 


the numerics. This phenomenon was overlooked in the recent paper 36 , since the focus was 
more on the physics about the S-dual (m = 0) point. 

The closing of the gap was expected since implementing the alternative quantization we 
allow the gauge held to be dynamical. In particular, this means that we let the external 
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Figure 10: We depict the (positive) real part of the lowest excitation mode in the D3-D5 
model with d = 10® and B = 3 ■ 10® as a function of ii. The dots are numer ically computed 
data points whereas the continuous curves stand for the analytic result of (6.15) for k = 0 
(lower curve) and k = 1 (upper curve). Notice that the mass gap for the zero sound vanishes 
precisely when ^ ^ for zero momentum. 


magnetic field to adjust its vacuum expectation value. The particular case where the gap 
closes corresponds to an anyonic fluid for which the effective magnetic field is vanishing 
and respectively the zero sound mode becomes gapless. Similar occurrence happened in the 
D3-D7’ model in the incompressible phase, where this mechanism was used to obtaining 
an anyonic superfluid, with the soft mode in the neutral sector becoming massless precisely 
when the effective magnetic held vanished 
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6.2 Conductivities 


Let us obtain the conductivity of the anyonic huid following the approach of 35 . As shown 


in that paper, the current-current correlator can be parametrized in terms of three functions 
Cl, Ct, and W, which transform under S-duality as: 

Cl Ct 


Cl = 


(27r)2 [ClCt +fF2] ’ 
W* = 


_ 

— 


(27r)2 [ClCt +fF2] ’ 


W 


(27r)2 [Cl Ct + ‘ 

Under the T transformation only the parity violation function W transforms as: 


(6.16) 


W* = IF + — 

27r 


(6.17) 
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The current-current correlator at zero momentum is parametrized in terms = 

Ct{oo,0) and lV(a;,0) as: 


{Ji{u,0) Jj{uj,0)) = iu [CL{u:,Q)5ij + , 

where i,j are spatial indices. The AC conductivities are dehned as: 


= — (Ji(a;,0) J,(a;,0)) , 

lOJ 


(6.18) 


(6.19) 


or, equivalently 

(Tl{u:) = Cl{u:,Q) , o-h(w) = W{u,0) . (6.20) 

In our case we can obtain the value of ctl at a; = 0 (the DC conductivity) by looking at the 
transverse correlator (C.41) at zero momentum. We get that an = 0 and that 

(Tl = , ((u ^ 0) , (6.21) 


which, after taking (C.42) into account, becomes: 




P+A-7 

^ M \J 1 + (P , (w —)■ 0) . 


( 6 . 22 ) 


Starting with a theory with W = 0 and performing the transformation we end up with 
a theory with a Hall conductivity an = K/{2 ti). With a subsequent S transformation the 
corresponding conductivities can be found from (6.16) with Ct(ci;,0 ) = = aiioj)-. 

1 / n 1 CrH(w) 


al{uj) = 


(27r) 


a?,(a;) + ( f 


2 ’ 


a*H{u) = 


+ {0' 


(6.23) 


For the AC conductivities there are no analytic results available in any model. Had there 
been for example results for the normal quantization, one could just transfer them to any 
quantization using the formulas (6.23). Given this shortcoming, we do not wish to make a 
thorough scan over all the available parameters in the model, but are elated to highlight just 
one interesting effect upon changing m. By numerically extracting the AC conductivities, 
as in 1^, we hnd that increasing m, the Drude-like peak will begin to transform to make a 


secondary peak at higher a), similar to what is found in heavy fermion systems; see Fig. n 
for how the longitudinal conductivities vary. 

The DC conductivities can be obtained by taking a; = 0 in the previous expression, as 
well as the value (6.22) for the DC conductivity ctl of the normally quantized system: 

01 + P 


(27r)^ r 


p+A-7 

2 - 2 ^ 


1 -i- ^2 -|- 


where, as in 35 , we identify m as: 


(2vr) r 


K 


p+A-7 

^ = 


m 


1 + d? + 0 


m = 


2nN ’ 


(6.24) 


(6.25) 


and m is related to m as in (6.9). 
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Figure 11: We are plotting the optical conductivities versus Cj for the D2-D6 model at hxed 

. ^ P+A-7 

d = 10 and i? = 0. In the left panel we are plotting the real part of (27r)^r^ ^ in 

the right panel the imaginary part. The different curves correspond to different quantization 
schemes m = 0,1, 2, 5,10. The bigger amplitudes correspond to smaller m’s. 

6.3 Diffusion constant 


Let us consider the current-current correlator with non-zero momentum k following the 
approach of [^. In this case we have to distinguish between longitudinal and transverse 
functions Cl and Ct- The longitudinal function Cl and its S-dual can be parametrized 
in terms of the DC conductivities ctl and and diffusion constants D and D* as: 


Cr = 


kaL 


iCj — Dk"^ 


Cl = 


k a*T 


iu — D* k? 


Notice that Cl and Cl are related as in (6.16), which we now write as: 

1 

= Ct + ^ . 


(27r)2 Cl 


Cr 


Moreover, the DC conductivities aL and al are related as: 

1 _ 

{2nYal ^ aL ' 


(6.26) 


(6.27) 


(6.28) 


Let us now obtain the transverse correlation function Ct in terms of Cl and Cl from (6.27). 
By using the parametrization (6.26) and the S-dual relation (6.28) we get that: 

ctl 


Ct= {toj 
k 


where H is given by: 


hP) , 

(6.29) 


(6.30) 


a 


Inverting this last relation we get D* in terms of H and D\ 

H + ^D 


D* = 


1 + 


n/2 


(6.31) 
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If parity is conserved in the original theory and, thus, W = 0, then H = D* and we can 
write Ct as: 

Ct = ^ {icu - D* k^) . (6.32) 

k 

We will use this relation to obtain D* from the transverse correlator calculated in section 0 
Indeed, when W = 0 the non-vanishing different components of (when k^ = {uj, k, 0) 
are related to and Ct as: 

7 2 2 

{Jt{k) u-k)) = - CT{k ), {Uk) u-k)) = - \ CT{k ), 

iUk) Jy{-k)) = ^L= CL{k) , iUk) Jyi-k)) = V-OJ^ + k^Crik) . 

y—oj^ + k^ 


(6.33) 


In the hydrodynamic regime in which k ~ 0(e) and u ~ at leading order, we get: 


{Uk)jy{-k))^kCT 


(6.34) 


Therefore, the expression of Ct can be found from the correlator (C.41). By comparing this 
equation with (6.32) and using (6.22), we hnd 

f, 2 1 


D* = 


1 p^l6 — p 16 — p 1 -2 

TTTJ ^ 2(6-p)-A 


(6.35) 


We can now use this value of D* to apply the argument of 
constant after a S transformation, which is given by: 


35 


D* + 


Dm = 


K 


D 


1 + 


K 

2-71 (TL 


and hnd the diffusion 


(6.36) 


which, with the identihcation (6.25) and the expression of ai written in (6.22), becomes: 

{l + d^)D* + m^£) 


Dm = 


1 -|- -|- lii^ 


(6.37) 


By using the values of D (eq. (4.24)) and D* (eq. (6.35)) we get: 


Dm = 


2 vTT^ 

1 + d'^ + m^ 


2(6-p)-A V2’ A 2’ A 2’ ^ 


m 


A-2 


3 1 13 


2’2 A’2 A’ 


-d^ 


(6.38) 
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Figure 12: The diffusion constants in the D2-D6 model with different quantization conditions. 
The red points represent numerical data. The different curves are the analytic results. The 
monotonically increasing solid blue curve corresponds to normal quantization, i.e., n = 0 
(4.24), the others to alternative quantization (6.38); the middle two dashed black curves 
correspond to m = 1000 , 100 , the dotted blue curve (visible as the middle on the left) to 
m = 1, and the monotonically decreasin g soli d blue curve to m = 0. Notice that all JDm 
(iff m < oo) asymptote to S-dual result (6.35) as d —)■ oo and at large temperatures they 
asymptote to constants as given in (6.39). 


From this formula we can get the limiting value of when d —0: 

2 \-7+p 


lim Dm = 

d^O 


A-2 (A-2)(12-2p-A) ' 


(6.39) 


The numerical values of Dm match extremely well with the analytic results (6.38) as repre¬ 
sented in Fig. f^for the D2-D6 model. 


7 Discussion: scaling behavior 


The energy scale of the bulk theory is given, in terms of the radial coordinate r, as 


c ~ r 2 


(7.1) 


The coefficient in (7.1) contains y/gvM, where qym is the Yang-Mills coupling of the bulk 
theory, which is dimensionful when p 7 ^ 3. This distance/energy relation was found long time 
ago in [^. It can also be derived by noticing that the Dp-brane background is conformally 
AdS. This means that one can perform a Weyl transformation to a dual frame: 


ds 


20 


dual 


= ds 


(7.2) 
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followed by a change in the radial coordinate of the form u ^ r~, after which the metric 
is that of AdSp +2 x The relation (7.1) follows by identifying the new radial coordinate 
u with the energy scale. We can also extract (7.1) by looking at the relation between the 
horizon radins Vh and T (2.18). 


Let ns now consider rescalings of S of the form: 


(7.3) 

In terms of p these rescalings are eqnivalent to: 

Ap=-^. (7.4) 

5 — p 

We will say that p transforms with scaling dimension Ap. 

Let ns obtain the behavior of the density d and magnetic held B by imposing that the 
three terms in the combination + p^+P“’^ appearing in the eqnation of motion for 

At (and also for the hnctnation helds), scale in the same way. We write 

d^A^^d, B^A^^B. (7.5) 


We hnd the following valnes for the scaling dimensions A^ and A^: 


Ad 


A 

5-p ’ 


As 


7 — p 
5 — p 


(7.6) 


Therefore, A determines the scaling dimension Ad of the charge density]^ while As depends 
only on p. In the conformal case p = 3 we have A = 2n and, thus. Ad = n and As = 2. These 
values are just the canonical dimensions of these helds in a (n + l)-dimensional QFT. In a 
non-conformal background Ad and As diher from the canonical dimensions. For example, 
for the supersymmetric intersections of Table we get: 

Af ^ (n + 3 - p) . (7.7) 

5 — p 5 — p 


Curiously, when p 7 ^ 3, = n only for (2 + l)-dimensional intersections (7e., with 

n = 2 for any p 7 ^ 3, corresponding to D2-D6 and D4-D4 in Table [^. 

The scaling dimensions of the diherent quantities determine the way in which they must 
be rescaled with the temperature to make them temperature independent. Indeed, given 
a quantity O, its rescaled quantity O is dehned in such a way that O is invariant under 
rescalings. As As = 1, it is clear that the relation between O and O must be of the type: 


O 


O 

T^o ’ 


(7.8) 


=’This was expected as we related A to the polytropic index of the equation of state, see (2.29). 
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with Aq being the scaling dimension of O. Notice that this agrees with o ur de hnition of u 
and k in (4.21) (since = Ah = 1) and also with the dehnition of d in (4.20). Recall the 
expression of B written in (5.8): 


B = 


B 


V 47r 


B 


7-p 


(7.9) 


A posteriori, we notice that the rescaling in (5.8) was chosen to agree with the general 


dehnition (7.8) when is given by (7.6) 


The behavior under rescalings can also be used to determine the general formula of the 
invariant quantities normalized by the density. These quantities will be denoted by a tilde 
and are generally given by: 

o = -T . (7.10) 

Let us hrst apply this dehnition to the temperature T, frequency u, and momentum k. We 
have: 


T = 


T 


d^d d A 


U} = 


Ul U} 

~A~ ~ ,5-p ) 

d^d d A 


k = 


k 


k 


J_ 5 ^ • 

d^d d A 


(7.11) 


For p = 3 the exponent of the density in these dehnition is 1/n, in agreement with the 

(7.12) 


dehnitions used in 22 . In the case of the magnetic held we have: 

B 


B = 


d A 


In the conformal case p = 3 the power of d is 2/n, again in agreement with (22). 


8 Conclusions and outlook 


In this paper we studied the collective excitations of holographic matter engineered as in¬ 
tersection of two stacks of D-branes. One type of branes (the color branes) were substituted 
by the geometry they generated, while the havor branes were considered as probes and their 
dynamics were governed by the DBI action. We analyzed these systems at high baryonic 
density, both at zero and non-zero temperature, and we also studied the inhuence of the 
magnetic held. 

We determined (both analytically and numerically) the dispersion relation of the holo¬ 
graphic zero sound at low temperature, as well as the dihusion constant at higher temper¬ 
ature. We also studied numerically the crossover between the collisionless regime at low 
temperature and the hydrodynamic regime at higher temperature. When the intersection 
is (2 -|- l)-dimensional, one can further study the anyonic degrees of freedom in the system 
by performing an alternative quantization. We implemented this procedure for a general 
(2 -|- l)-dimensional brane intersection and computed the corresponding anyonic correlators. 
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Our results apply to a large number of brane intersections, characterized by the three 
numbers p, q, and n. However, we found that they only depend on p (the dimensionality of the 
background branes) and the index A dehned in (1.1). Indeed, this (p. A) universality showed 
up in the different rescalings we performed and in the relations that the rescaled quantities 
satished. Thus, our analysis unihed and extended previous results in the literature. 


We are continuing our efforts in generalizing the results presented here to the case, where 
we allow a non-zero mass for the fundamentals. This is a much more involved study, since 
then also the fluctuations of the scalar held has to be taken care of. However, our initial 
exploration suggest that most of the results in the present paper can be made available also in 
the presence of the mass. Actually, for the supersymmetric intersections at zero temperature 
and nonvanishing chemical potential, one can choose a system of coordinates such that the 
embedding function of the probe is also a cyclic variable [^. This has allowed to compute 
analytically the dispersion relation of the zero sound for massive quarks in the case of the 
D3-D7 10 19 , D3-D5, and D3-D3 intersections j^. Interestingly, the speed of zero sound 
for these systems vanishes when the quark mass equals the chemical potential, signaling a 
phase transition with violation of hyperscaling, whose critical exponents have been evaluated 
We intend to analyze this phenomenon for arbitrary intersections in the near future. 


m 


60 


Another interesting avenue to pursue would be to allow an internal flux in the worldvol- 
ume of the probe D-brane. In the Dp-D(p-|- 2) systems, turning on the appropriate internal 
flux induces a non-trivial prohle of some scalar. This corresponds, in the held theory dual, 
to moving to the Higgs branch of the theory (see [^). The analysis of the collective exci¬ 
tations of these intersections is an interesting open problem which we will try to address in 
the future. 
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A Fluctuation equations of motion 


Let us write the equations of motion (3.15) with B held. The equation of motion for at 
becomes: 


dn 


Idttl H j Qxx fp 


1 

+ “ 
U 


9rr H 


9x 


fp \9lx + B^ 


dx {dx Oit dt Ux) 


+do 


dB 


9lx + B^ 


dx ay 0 


(A.l) 


Let us write this equation in momentum space. By using (3.21) we can write it in terms of 
the electric held E. The result is: 


E" + a, log 


I 9xx fp ^H -|- d 


^ V9rr 9lx + B^ UJ^- /c2 ^2 
= i B d 


E' + 


9r 


\9tt\ fp 


{u;^-k^u^)E 


gfx + B^ 


dn 


\/\gtt\ 9xx fp y/H + I 5'xx + 


“y • 


(A. 2 ) 


For a Dp-brane background the metric elements and the blackening factor are given in (2.16) 
and (2.17) while the function H has been written in (2.19). Moreover, when i? 7 ^ 0 we have 

(A.3) 


'Z 7 2 2 

OJ — k U = 


{up — fp k'^) + up d^ 


pX _|_ pX+p-7 Q2 _|_ (P 

Plugging these values in the previous equation, we get: 

E~P {p^ + p^+P-7 i^2 _|_ I 


E" + dp\og 


P 


p7-p + 52 (^2 _ p ^2) pX + ^2 52 pX+p-7 + ^2 p 
1 (a;2 - fp e) p^ + a;2 B^ p^+P-'^ + u^d^ 


E' 


p7-p J2 pX _j_ pX+p-7 ]^2 _j_ (P 

{7-p)iBd {UJ^- fp e) p^ + UJ^ 52 p^+P-7 + a;2 d^ 


pfp + B^) 


{p^ + p^^P~'^ i?2 + d'^) 2 


"y • 


(A.4) 


The equation for ay can be written as: 


a'' + dp log 


9xx y/lfi'ul r pH + <72 

h 


9xx + B'^ ^ /2 \gtt 


9r 


(a ;2 — fc 2 


p 

2 I d2 


= -idB 


9 XX I f p \/H -\- d 


52 


E 


(A.5) 
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For the Dp-brane background this equation becomes: 


+ dp log 


P 


7 -p 




^pX + pX+P-7B^ + d^fp 


1 (a;2 - fp k^) 52 ^A+p -7 ^ ^2 p 2 

^p7-p/2 pA + ^A+p-7 52 + ^2 

{7 — p)iBd E 

Pfp {p'^-P + 52) ^pA + pA+p -7 52 + ^2 


(A.6) 


B Indicial equation 


Let us consider a differential equation for the function E = E{p) of the type: 

1 


E" + 


(- + Cl) E' + ( 

\p-rh / V 


A _|_ *^2 


(p - PhY p-Ph 


E = 0 , 


(B.l) 


where r/j, A, ci, and C 2 are constants. We want to hnd a solution in Frobenius series for p 
close to Th, of the type: 


E{p) = Enh {p - PhT 1 + Y{p-Ph) + • • • 


(B.2) 


By substituting this expansion in (B.l) and comparing the different powers of p — Vh, we get 


that a and fd must satisfy the equations: 


+ A = 0 , 

“h 2o( -|- -|-1) "1“ OLC\ -f- C 2 — 0 . 

(B.3) 

By choosing infalling boundary conditions, we get that a (when A > 0) must be: 



a = —i \/A . 

(B.4) 

Moreover, jd is given by: 

aci + C 2 
^ ~ 1 + 2a 

(B.5) 


C Transverse correlators 


We now study the equation of motion (4.2) for ay in the low frequency regime in which 
k ~ 0(e) and u ~ 0{e‘^). First, we expand (4.2) near the horizon p = Vh- The coefficient 


of the term without derivatives is just the same as in (4.1) and can be expanded as in (4.4) 
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and (4.5). Moreover, the coefficient multiplying a' in (4.2) can be represented near p = 


as: 


where di is given by: 


dp log 

^/p^T^ fp 

1 

— -|- di -|- ... , 

(C.l) 



p-Th 



V + X — 

8 A d^ 

(C.2) 

di — 

2rh d? + r^ 

in Frobenius series 

around p = Th. 


ay{p) 

= {p-rhT 

{1 +P{p-rh) + ...) . 

(C.3) 


From the equations of appendix one can show that the exponent a is just the 


same as m 


(4.6). The coefficient [3 is given by (B.5) with ci changed by di (and given by (C.2)). Since 
a, C 2 ~ 0{e‘^) and di ~ 0{1), we have that, at order e^, [3 can be written as: 


(3 ^ -{adi + C 2 ) . 


(C.4) 


Plugging the values of a, di, and C 2 (written in ( |4.6 ), (C.2), and (4.5) respectively), we hnd 
that (3 is given by: 


P = 


2(7-p) 


7 -p 


(p - 8) d^ + (p + A - 8) 

d^ + rl 


OJ + 


r 


,p+A-6 

h 


{7-p) {d^ + r^) 


e . 


(C.5) 


We now take the near-horizon and low-frequency limits in the opposite order. First of all, 
let us write the equation of motion of ay in the form: 


G' 

Oy + ~Q % + Q tty = 0 , 


where G = G{p) is dehned as: 


G{p) = Vp^ + d2/p(p) 


(C.6) 


(C.7) 


The function Q in (C.6) can be read from (4.2). Notice that, at order e^, this function is 
simply: 


Qip) 


„\+p —7 


{p^ + d2) fp{p) 


(C.8) 


We want to match the near-horizon expansion (C.3). Therefore, it is better to redehne ay in 
the form: 


The resulting equation for Uy is: 


G^ 

~G 


ay(p) = F{p) ^yi.p) , 

(C.9) 

= Th and F{p) is given by: 


F{p) = (P- rnT ■ 

(C.IO) 

2 — j a'y + e'^ {P + Q) ay = 0 , 

(C.ll) 
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where we have explicitly introduced the powers of e to keep track of the low frequency 
expansion and we have dehned the new function P{p) as: 



77// yO/ JPf 

(C.12) 

We will solve (C.ll) order by order in a series expansion in e of the form: 



- 0^0 ^ “!“••• • 

(C.13) 

The equation for ckq is: 

G' 

“0 + “0 = 0 ) 

(C.14) 

whose hrst integration yields: 


/ <-0 

(C.15) 


where cq is a constant. When cq 7^ 0 the function a'^ blows up at p = Vh- Thus, we choose 
Co = 0 and therefore ao = constant. Without loss of generality we can take 


The equation for ai is 


ao — 1 


a'[ + — a'l — —P — Q 

Cj 


Let us solve this equation by variation of constants. We put: 

A(p) A(p) 


a^ = 


\/p^ + (P fp{p) G{p) ’ 


(C.16) 

(C.17) 

(C.18) 


where A(p) is a function to be determined. By direct substitution into (C.17) we get that 
A(p) must satisfy: 

A' = -G(P + Q). (C.19) 

To integrate this equation we notice that the hrst term on the right-hand-side is, actually, a 
total derivative. Indeed, at 0{e^) we have that 


(logF) = — - - « — , 


(C.20) 


since F'/F ~ when F is the function (C.IO). Using this result in the dehnition of P 

we conclude that: 


GP = 


G(logF)' +0{e^). 


Therefore, 


A(p) = -G(logF)' - c 


G{p)Q{p) dp , 


'rn 


where c is a constant to be determined. Let us next dehne the integral X(p) as: 

k^Kp) = - f G(p)Q(p) dp, 

J rh 


(C.21) 

(C.22) 

(C.23) 
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or, more explicitly: 


Then, ol^ can be written as: 


i(p) 



pA+p—7 


dp . 


(C.24) 


tti = 


+ 


a 


\/p^ + d? fp P-Th 


+ 


Ap) 

fp 


(C.25) 


We will fix the constant c by imposing that a'^ be regular at p = r^. This is equivalent to 
require that the term with c cancels the term with l/(p — rh)- Notice that, by construction, 
X(p) vanishes linearly at the horizon and the last term in ( |C.25 ) is therefore regular at the 
horizon. Taking into account the expansion of l/G around p = r/^: 


rh 


{p — 8)d‘^ + (p + A — 8) r 


AP^ + d‘^ fp {7 - p) ^/r^A~cP P - rh 2(7 - p)(r^ + 2 


we get: 


c = 


7 — p 
rh 


A + d? a = i 


Ar^ + d‘^ 


7 -p 
, 2 


UJ . 


(C.26) 


(C.27) 


Thus, a' = a'l + 0{e^) and is given by: 


A = 


(7 - p) 


5 -p 


7-p y/r^ + d2 


rh A P^ + fp 


CO + 


Ap) 


Ap^ + dP fi 




(C.28) 


In order to match this solution with ( |C.3 ), let us expand X(p) near p = rh- We easily get: 

{p - Th) + 0{{p - Thf) - (C.29) 


Ap) = 


r^+P-7 
' h 

A A + 


From this result and the expansion of 1/G in (C.26), it is easy to check that, indeed, 

= rh) = A (C.30) 


where jd is given in (C.5). Moreover, since ~ 1 + ai and F 1 + a log(p — Vh), we 
conclude that, at (!7(e^): 


tty ~ 


1 + 0:1 + Olog(p - Th) 


Therefore, it follows that: 


Explicitly, 


f f I 

dy Ol + 


+ , = 


+ oA) = o' + + oA) 

p-rh ^ p-rh 


1 . AA + d? n-( \ 1-2 

-- 777 -. ^ — u-I{p)k 

G{p) 


(C.31) 


(C.32) 


(C.33) 
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where G{p) has been dehned in (C.7). 


Let us now extract the {Jy Jy) correlator from the previous results. Recall that the term 
depending on ay of the Lagrangian density is of the form: 


£K) = ■ 


where is given by: 


j, ^ gyy gpp ^ G{p) , 

VH + d^ ^ ’ 

with J\f being a normalization constant. The on-shell action of ay is: 


*S^on—shell(^^) j d X 


p^oo 


(C.34) 


(C.35) 


(C.36) 


With the normalization condition we are using (ao = 1), the two-point function of Jy is 
given by the standard AdS/CFT prescription: 


{Jy{p) Jy{-p)) = 

From the explicit expressions of and a^, we get: 

^rl + 


p^OO 


(C.37) 


Tdy = N 


T-V 
. 2 


oj — I{p) 


(C.38) 


In order to get the value of the integral X at the boundary, let us rewrite its expression as: 


rp’’ ° Ap) ( 

i - 


p—6 

A 


dC, 


(C.39) 


C + d^ 


with d given by (4.20). Then, it is easy to prove that, when 2(6 — p) > X the p —)■ oo limit 

l+p-6 6-p 16-p 1 


of X converges and is given by: 

2 


X(p —)■ cx)) = 


F p, 


2 (6 -p)-A ^ 

It follows that the {Jy Jy) correlator takes the form: 


+ -■ 

2’ A 2’ A 2’ 


{Jy{p) Jy{-p)) = M 
where the coefficients F(^ and F^ are: 




(C.40) 


(C.41) 


A+p-7 

F^ = \ll + d? 


F. = 


2(6-p) - A 




2’ A 2’ A 


(C.42) 
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When p = 3 and A = 4 the result written above coincides with the one in 35 . 


Let us rewrite the correlator in terms of the rescaled frequency and momentum uj and k 

(C.43) 


defined in (4.21): 


{Jyip) Jy{-p)) = N 


T^iuj + Tkk 


where N is related to N as: 

and the rescaled coefficients and are: 


M = M r'^ 


-1 


r. 

tk 


= \li + (P 
2 


2 ( 6 -p)-A V 2 ’ A 


(I 6-p _ 1 6-p 1 

2 ’ A ^2’ 


(C.44) 


(C.45) 


D Wronskian method 


In this appendix we solve the inhomogeneous equation (5.25) by the Wronskian method. We 
start by defining a new function y as: 


y = P ^ y{p) , 


(D.l) 


and a new independent variable x as: 


2cu p—5 

X = - p 2 

5 — p 


(D.2) 


In what follows we consider 3^ as a function of x. After this change of variables, eq. (5.25) 
becomes: 

d^y 1 dy 


dx"^ 


X dx 


where f{x) is given by: 


and Up and dp are defined as: 


Up = 


5 — p 



(D.3) 


(D.4) 

0--(7-p)(=7)*^-^. 

(D.5) 


Eq. (D.3) is a linear inhomogeneous equation, whose solutions can be obtained from those 
of the homogeneous equation. Let and y 2 {x) be two independent solutions of (D.3) 

with f{x) = 0. Then, the solution of (D.3) for f{x) ^ 0 can be written as: 


y{x) = h{x)yi{x) + i2{x)y2{x ), 


(D,6) 
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where Ii{x) and l 2 {x) are the following indehnite integrals: 


lAx) = - 


y2{x)f{x) 


dx , 


h{x) = 


and hh(3^i, 3 ^ 2 ) is the Wronskian function of and 3 ^ 2 : 

hk( 3 ^i, 3 ^ 2 ) = 3^13^2' - 3^23^1' . 


yi{x)f{x) 

w{y^,y2) 


dx , 


(D,7) 


(D,8) 


The homogeneous version of (D.3) is just the Bessel equation. Therefore, we can take the 
Hankel functions of index Vp as the two independent solutions 3^i and 3 ^ 2 : 


= (i = l,2). 

The Wronskian of two Hankel functions is rather simple, namely: 

Ai 

w(h<i\x),hS,>m) = —. 


TTX 


(D,9) 


(DIO) 


Therefore, Ii(x) and hix) are given by: 

r / \ • 

hix) = -l—r- 


.'Xdr, 


hix) = i- 


Taking into account that (for z/ + /i 7 ^ 1): 


j (a^) hII^ ( 2 :) dx 

x'^'^p+^ hII\x) Hl^J{x)dx . 


(D.ll) 


xi^+^+^ Hjy\x) Hjf\x) dx = 


X' 




j 2(/i + 1/ + 1) L ^ 

we get that hix) and hi^x) are given by: 


H'fHx) HfAix) 


(D.12) 


/i(a;) = -i — 


.Tidp 5 — p 2(6-P) 


7 — p 


X 


H!:^Hx) Hf'Hx) + H^^Ux) H<li,(x) 


r(l) 


l2(x) = i 


It then follows that: 


. 71 dp 5 — P 2(6-p) 


7 — p 


■ X 5-P 


H^Hx) Hll\x) + HilA(x) XJAdA ] ■ (D.13) 


( 1 ) 

'p +1 '' 


y{x) = i 


. ndp 5 — p 2(6-p) 


7 — p 


X 5-p 


<b(p) 1 <h(p) ■ (D-14) 


By using the following property of the Hankel functions: 


- h !.^’( x ) h 7+,( x ) = -y, 

X 


( 2 ) 

'p- 


Ai 


we arrive at: 


y{x) = 


{5-p)x]i4 Ml) 


2u 




5-p 


(D.15) 

(D.16) 


which matches with the expression for y{p) in (5.26). 


57 




























References 


[1] L. D. Landau, “The theory of a Fermi liquid”, Zh. Eksp. Teor. Fiz. 30, 1058 (1956) 
[Soviet Phys. JETP 3, 920 (1957)]. 

[2] See, for example: E. M. Lifshitz, L. P. Pitaevskii, “Statistical Physics”, Part 2, Perg- 
amon Press, Oxford 1980; D.Pines and P. Nozieres, “The theory of quantum liquids”, 
Benjamin, New York 1966. 

[3] For reviews see: J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal and 
U. A. Wiedemann, “Gauge/String Duality, Hot QCD and Heavy Ion Collisions,” 
arXiv: 1101.0618 [hep-thj; J. McGreevy, “Holographic duality with a view toward many- 
body physics,” Adv. High Energy Phys. 2010, 723105 (2010) |arXiv:0909.0518 [hep-thj]; 
A. V. Ramallo, “Introduction to the AdS/CFT correspondence,” Springer Proc. Phys. 
161 (2015) 411 |arXiv:1310.4319 [hep-th]]. 

[4] A. Karch and E. Katz, “Adding flavor to AdS / CFT,” JHEP 0206 (2002) 043 |hep- 
th/0205236|. 

[5] S. Kobayashi, D. Mateos, S. Matsuura, R. C. Myers and R. M. Thomson, “Holographic 
phase transitions at finite baryon density,” JHEP 0702 (2007) 016 |hep-th/0611099|. 

[6] A. Karch, D. T. Son and A. O. Starinets, “Zero Sound from Holography,” 
arXiv:0806.3796 [hep-th]. 

[7] A. Karch, D. T. Son and A. O. Starinets, “Holographic Quantum Liquid,” Phys. Rev. 
Lett. 102 (2009) 051602. 

[8] O. Bergman, N. Jokela, G. Lifschytz and M. Lippert, “Striped instability of a holo¬ 
graphic Fermi-like liquid,” JHEP 1110 (2011) 034 |arXiv:1106.3883 [hep-th]]. 

[9] N. Jokela, G. Lifschytz and M. Lippert, “Magnetic effects in a holographic Fermi-like 
liquid,” JHEP 1205 (2012) 105 |arXiv:1204.3914 [hep-th]]. 

[10] M. Kulaxizi and A. Parnachev, “Comments on Fermi Liquid from Holography,” Phys. 
Rev. D 78 (2008) 086004 |arXiv:0808.3953 [hep-th]]. 

[11] K. Y. Kim and 1. Zahed, “Baryonic Response of Dense Holographic QCD,” JHEP 0812 
(2008) 075 |arXiv:0811.0184 [hep-th]]. 

[12] M. Kulaxizi and A. Parnachev, “Holographic Responses of Fermion Matter,” Nucl. 
Phys. B 815, 125 (2009) |arXiv:0811.2262 [hep-th]]. 

[13] L. Y. Hung and A. Sinha, “Holographic quantum liquids in Ij-l dimensions,” JHEP 
1001 (2010) 114 |arXiv:0909.3526 [hep-th]]. 

[14] M. Edalati, J. 1. Jottar and R. G. Leigh, “Holography and the sound of criticality,” 
JHEP 1010 (2010) 058 |arXiv: 1005.4075 [hep-th]]. 


58 





[15] B. H. Lee and D. W. Pang, “Notes on Properties of Holographic Strange Metals,” Phys. 
Rev. D 82 (2010) 104011 |arXiv:1006.4915 [hep-th]]. 

[16] C. Hoyos-Badajoz, A. O’Bannon and J. M. S. Wn, “Zero Sonnd in Strange Metallic 
Holography,” JHEP 1009 (2010) 086 |arXiv:1007.0590 [hep-th]]. 

[17] B. H. Lee, D. W. Pang and C. Park, “Zero Sonnd in Effective Holographic Theories,” 
JHEP 1011 (2010) 120 |arXiv:1009.3966 [hep-th]]. 

[18] M. Ammon, J. Erdmenger, S. Lin, S. Mnller, A. O’Bannon, J. P. Shock, J. Erdmenger 
and S. Lin et al, “On Stability and Transport of Cold Holographic Matter,” JHEP 
1109 (2011) 030 |arXiv:1108.1798 [hep-th]]. 

[19] R. A. Davison and A. O. Starinets, “Holographic zero sonnd at hnite temperature,” 
Phys. Rev. D 85 (2012) 026004 |arXiv: 1109.6343 [hep-th]]. 

[20] M. Goykhman, A. Parnachev and J. Zaanen, “Fluctuations in hnite density holographic 
quantum liquids,” JHEP 1210 (2012) 045 |arXiv:1204.6232 [hep-th]]. 

[21] A. Gorsky and A. V. Zayakin, “Anomalous Zero Sound,” JHEP 1302 (2013) 124 
|arXiv: 1206.4725 [hep-th]]. 

[22] D. K. Brattan, R. A. Davison, S. A. Gentle and A. O’Bannon, “Collective Excita¬ 
tions of Holographic Quantum Liquids in a Magnetic Field,” JHEP 1211 (2012) 084 
larXiv: 1209.0009 [hep-th]]. 

[23] N. Jokela, M. Jarvinen and M. Lippert, “Fluctuations and instabilities of a holographic 
metal,” JHEP 1302 (2013) 007 |arXiv:1211.1381 [hep-th]]. 

[24] R. A. Davison and A. Parnachev, “Hydrodynamics of cold holographic matter,” JHEP 
1306 (2013) 100 |arXiv:1303.6334 [hep-th]]. 

[25] D. W. Pang, “Probing holographic semilocal quantum liquids with D-branes,” Phys. 
Rev. D 88 (2013) 4, 046002 |arXiv:1306.3816 [hep-th]]. 

[26] P. Dey and S. Roy, “Zero sound in strange metals with hyperscaling violation from 
holography,” Phys. Rev. D 88 (2013) 046010 |arXiv:1307.0195 [hep-th]]. 

[27] M. Edalati and J. F. Pedraza, “Aspects of Current Correlators in Holographic Theories 
with Hyperscaling Violation,” Phys. Rev. D 88 (2013) 086004 |arXiv:1307.0808 [hep- 
th]]. 

[28] R. A. Davison, M. Goykhman and A. Parnachev, “AdS/CFT and Landau Fermi liq¬ 
uids,” JHEP 1407 (2014) 109 |arXiv: 1312.0463 [hep-th]]. 

[29] B. S. DiNunno, M. Ihl, N. Jokela and J. F. Pedraza, “Holographic zero sound at hnite 
temperature in the Sakai-Sugimoto model,” JHEP 1404 (2014) 149 |arXiv: 1403.1827 
[hep-th]]. 


59 



[30] R. E. Arias and I. S. Landea, “Hydrodynamic Modes of a holographic p— wave super¬ 
fluid,” JHEP 1411 (2014) 047 |arXiv:1409.6357 [hep-th]]. 

[31] K. Jensen, A. Karch and E. G. Thompson, “A Holographic Quantum Critical Point 
at Finite Magnetic Field and Finite Density,” JHEP 1005 (2010) 015 |arXiv:1002.2447 
[hep-th]]. 

[32] W. Kohn, “Cyclotron resonance and the Haas-van Alphen oscillations of an interacting 
electron gas”, Phys. Rev. 123 (1961) 1242. 

[33] N. Jokela, G. Lifschytz and M. Lippert, “Holographic anyonic superfluidity,” JHEP 
1310 (2013) 014 |arXiv:1307.6336 [hep-th]]. 

[34] N. Jokela, G. Lifschytz and M. Lippert, “Flowing holographic anyonic superfluid,” JHEP 
1410 (2014) 21 |arXiv:1407.3794 [hep-th]]. 

[35] D. K. Brattan and G. Lifschytz, “Holographic plasma and anyonic fluids,” JHEP 1402 
(2014) 090 |arXiv:1310.2610 [hep-th]]. 

[36] D. K. Brattan, “A strongly coupled anyon material,” arXiv: 1412.1489 [hep-th]. 

[37] A. Karch, M. Kulaxizi and A. Parnachev, “Notes on Properties of Holographic Matter,” 
JHEP 0911 (2009) 017 |arXiv:0908.3493 [hep-th]]. 

[38] S. A. Hartnoll, J. Polchinski, E. Silverstein and D. Tong, “Towards strange metallic 
holography,” JHEP 1004 (2010) 120 |arXiv:0912.1061 [hep-th]]. 

[39] F. Bigazzi, A. L. Cotrone and J. Tarrio, “Charged D3-D7 plasmas: novel solutions, 
extremality and stability issues,” JHEP 1307 (2013) 074 |arXiv: 1304.4802 [hep-th]]. 

[40] D. Arean and A. V. Ramallo, “Open string modes at brane intersections,” JHEP 0604 
(2006) 037 |hep-th/0602174|. 

[41] T. Sakai and S. Sugimoto, “Low energy hadron physics in holographic QCD,” Prog. 
Theor. Phys. 113 (2005) 843 |hep-th/0412141|. 

[42] O. Bergman, N. Jokela, G. Lifschytz and M. Lippert, “Quantum Hall Effect in a Holo¬ 
graphic Model,” JHEP 1010 (2010) 063 |arXiv:1003.4965 [hep-th]]. 

[43] R. C. Myers and M. C. Wapler, “Transport Properties of Holographic Defects,” JHEP 
0812 (2008) 115 |arXiv:0811.0480 [hep-th]]. 

[44] N. Jokela, M. Jarvinen and M. Lippert, “A holographic quantum Hall model at integer 
hlling,” JHEP 1105 (2011) 101 |arXiv: 1101.3329 [hep-th]]. 

[45] G. Itsios, N. Jokela and A. V. Ramallo, “Cold holographic matter in the Higgs branch,” 
Phys. Lett. B 747 (2015) 229 |arXiv:1505.02629 [hep-th]]. 

[46] P. M. Hohler and M. A. Stephanov, “Holography and the speed of sound at high tem¬ 
peratures,” Phys. Rev. D 80, 066002 (2009) |arXiv:0905.0900 [hep-th]]. 


60 



[47] A. Gherman, T. D. Cohen and A. Nellore, “A Bound on the speed of sound from 
holography,” Phys. Rev. D 80 (2009) 066003 |arXiv:0905.0903 [hep-th]]. 

[48] A. Gherman and A. Nellore, “Universal relations of transport coefficients from hologra¬ 
phy,” Phys. Rev. D 80 (2009) 066006 |arXiv:0905.2969 [hep-th]]. 

[49] I. Amado, M. Kaminski and K. Landsteiner, “Hydrodynamics of Holographic Super¬ 
conductors,” JHEP 0905 (2009) 021 |arXiv:0903.2209 [hep-th]]. 

[50] M. Kaminski, K. Landsteiner, J. Mas, J. P. Shock and J. Tarrio, “Holographic Operator 
Mixing and Quasinormal Modes on the Brane,” JHEP 1002 (2010) 021 |arXiv:0911.3610 
[hep-th]]. 

[51] S. A. Hartnoll, “Theory of universal incoherent metallic transport,” Nature Phys. 11 
(2015) 54 |arXiv:1405.3651 [cond-mat.str-el]]. 

[52] E. Witten, “SL(2,Z) action on three-dimensional conformal held theories with Abelian 
symmetry,” In *Shifman, M. (ed.) et ah: From helds to strings, vol. 2* 1173-1200 |hep- 
th/0307041|. 

[53] H. U. Yee, “A Note on AdS / GET dual of SL(2,Z) action on 3-D conformal held theories 
with U(l) symmetry,” Phys. Lett. B 598 (2004) 139 |hep-th/0402115|. 

[54] K. Goldstein, N. lizuka, S. Kachru, S. Prakash, S. P. Trivedi and A. Westphal, “Hologra¬ 
phy of Dyonic Dilaton Black Branes,” JHEP 1010 (2010) 027 |arXiv:1007.2490 [hep-th]]. 

[55] A. Bayntun, C. P. Burgess, B. P. Dolan and S. S. Lee, “AdS/QHE: Towards a Holo¬ 
graphic Description of Quantum Hall Experiments,” New J. Phys. 13 (2011) 035012 
|arXiv:1008.1917 [hep-th]]. 

[56] M. Fujita, M. Kaminski and A. Karch, “SL(2,Z) Action on AdS/BCFT and Hall Con¬ 
ductivities,” JHEP 1207 (2012) 150 |arXiv:1204.0012 [hep-th]]. 

[57] M. Lippert, R. Meyer and A. Taliotis, “A holographic model for the fractional quantum 
Hall ehect,” JHEP 1501 (2015) 023 jarXiv: 1409.1369 [hep-th]]. 

[58] A. W. Peet and J. Polchinski, “UV / IR relations in AdS dynamics,” Phys. Rev. D 59 
(1999) 065011 |hep-th/9809022]. 

[59] A. Karch and A. O’Bannon, “Holographic thermodynamics at hnite baryon density: 
Some exact results,” JHEP 0711 (2007) 074 |arXiv:0709.0570 [hep-th]]. 

[60] M. Ammon, M. Kaminski and A. Karch, “Hyperscaling-Violation on Probe D-Branes,” 
JHEP 1211 (2012) 028 |arXiv: 1207.1726 [hep-th]]. 

[61] D. Arean, A. V. Ramallo and D. Rodriguez-Gomez, “Mesons and Higgs branch in defect 
theories,” Phys. Lett. B 641 (2006) 393 |hep-th/0609010|; “Holographic flavor on the 
Higgs branch,” JHEP 0705 (2007) 044 |hep-th/0703094, 


61 





